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ABSTRACT 

The object of the research reported herein was to develop a 
general mathematical model and solution methodologies for 
analyzing the structural response of thin, metallic shell 
structures under large transient, cyclic, or static 
thermomechanical loads. Among the system responses associated 
with these loads and conditions are thermal buckling, creep 
buckling, and ratcheting. Thus geometric and material 
nonlinearities (of high order) can be anticipated and must be 
considered in developing the mathematical model. The methodology 
is demonstrated through different problems of extension, shear . 
and of planar curved beam. Moreover, importance of the inclusion 
of large strains is clearly demonstrated, through the chosen 
applications . 
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1 . INTRODUCTION 


1,1 Background 

The prediction of inelastic behavior of metallic materials at elevated 
temperatures has increased in importance in recent years. The operating 
conditions within the hot section of a rocket motor or a modern gas turbine 
engine present an extremely harsh thermo— mechanlal environment. Large 
thermal transients are induced each time the engine is started or shut down. 
Additional thermal transients from an elevated ambient occur, whenever the 
engine power level is adjusted to meet flight requirements. The structural 
elements employed to construct such hot sections, as well as any engine 
component located therein, must be capable of withstanding such extreme 
conditions. Failure of a component would, due to the critical nature of 
the hot section, lead to an immediate and catastrophic loss in power and 
thus cannot be tolerated. Consequently, assuring satisfactory long term 
performance for such components is a major concern for the designer . 

Traditionally, this requirement for long term durability has been a 
more significant concern for gas turbine engines rather than rocket motors. 
However, with the advent of reusable space vehicles, such as the Space 
Shuttle, the requirement to accurately predict future performance following 
repeated elevated temperature operations must now be extended to includ the 
more extreme rocket motor application. These operating blades to severe 
thermal transients that result in large inelastic strains, and several 
types of behavior must be considered. The elevated temperatures can lead 
to thermal buckling and, in addition, creep can be expected to occur. 
Thus, a combination of thermal-creep buckling behavior leading to large 
deflections can be anticipated. Because of the cyclic character of the 
mechanical and thermal loads, progressive growth or ratchetting effects 
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must also be considered. Thus, geometric and material nonlinearities (of 
high orders) can be anticipated and must be considered in the mathematical 
model . 

Consequently, the industry is concerned with the behavior of thin 
shell like structural elements subjected to severe time dependent 
thermo-mechanical loading. Such thin elements, including beams, rings, 
arches, plates and shells, are presenting generic types of components, 

/ 

which might be located within or adjacent to the hot section of a rocket or 
a gas turbine engine. 

The experience in the gas turbine engine industry indicates, however, 
that existing analytic tools are not sufficiently reliable to accomplish 
this task. State of the art methods for predicting hot section component 
behavior are generally not sufficiently accurate to perform extended use 
evaluations. 

Under this kind of severe loading conditons, the structural behavior 
is highly nonlinear due to the combined action of geometrical and physical 
nonlinearities. On one side, finite deformation in a stressed structure 
introduces nonlinear geometric effects* On the other side, physical 
nonlinearities arise even in small strain regimes, whereby inelastic 
phenomena play a particularly important role. From a theoretical 
standpoint, nonlinear constitutive equations should be applied only in 
connection with nonlinear transformation measures (implying both 
deformation and rotations). However, in almost all of the works in this 
area, the two identified sources of nonlinearities are always separated. 
This separation yields, at one end of the spectrum, problem of large 
response, while at the other end, problems of viscous and/or non- isothermal 
behavior in the presence of small strain. 
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Because of the nature of the causes, special care is needed in the 
selection or development of a constitutive law that Includes tme and 
temperature effects. Although there exists a sizeable body of literature 
on phenomenological consitutive equations for the rate - and 
temper ature - dependent plastic - def ormation of metallic materials, to date 
rational and thermodynamically consistent elastic - thermovlscoplastic 
constitutive relations capable of incorporating the effects of large 
strains and rotations have not been demonstrated. 

Constitutive models for small strain in engineering literature may 
generally be grouped into three categories: classical plasticity, 

nonlinear viscoelasticity, and theories based on microstructural phenomena. 
Each group can be further separated into "unified" and "uncoupled" 
theories, where the two differ in their approach to the treatment of rate - 
independent and rate-dependent inelastic deformation. The uncoupled 
theories decompose the inelastic strain rate into a time-independent 
plastic strain rate and a timedependent creep rate with independent 
constitutive relations describing plastic and creep behavior. Such 
uncoupling of the strain components provides for simpler theories to be 
developed but precludes any creep-plasticity interaction. Recognizing that 
cyclic plasticity .creep and recovery are not independent phenomena but 
rather are very interdependent, a number of "unified" models for inherently 
time- dependent nonelastic deformation have been developed recently. 

Classical incremental plasticity theories are macrophenomenological 
because they base the derivation of state variables purely on experimental 
results without direct reference to the micro3tructure of the material. 
Most incremental plasticity theories have four major components: (1) a 

stress-elastic strain relations, (2) a yield function describing the onset 
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of plastic deformation, (3) a hardening rule which prescribes the 
strain-hardening of the material and the modification of the yield surface 
during plastic flow, and ( 4 ) a flow rule which defines the components of 
strain that are plastic or nonrecoverable. Research in this area is 
voluminous. For example, Zienkiewicz and Cormeau^ developed a rate 
dependent unified theory which allows for nonassoclati ve plasticity and 
strain softening, but does not model the Bauschinger effect or temperature 
dependence. Extensions of classical plasticity to model both rate and 
temperature effects were presented recently by Allen and Haisler 2 , Haisler 
and Cronenworth3, and Yamada and Sakurai^. 

In the nonlinear viscoelastic approach, the constitutive relation is 
expressed as a single integral or convoluted form. This type of 
constitutive model employs the thermodynamic laws along with physical 
constraints to complete the formulation. A detailed review of several 
existing theories is presented by Walker^. Walker's® theory is based on a 
unified viscoplastic integral developed by modifying the constitutive 
relations for a linear three paramter viscoelastic solid. The theory 
contains clearly defined material parameters, a rate dependent equilibrium 
stress, and a proposed multiaxial model. An important shortcoming of 
Walker's theory is its failure to model transient temperature conditions. 
Many other nonlinear viscoelastic theories have been proposed including 
those by Cernocky and Krempl®, Valanis? and Chabache®. 

The microphenomenological theories attempt to represent the response 
of polycrystalline materials in terms of various micromechanisms of 
deformation and failure. Various dislocation theories have been developed 
to predict plastic deformation in terms of dislocation interaction, slip, 
glide, density, etc. Most of the material models developed, to date, 
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depend primarily on the number of state variables used and their growth or 
evolutionary laws. Many of the recent "unified" ml crophenomenologi cal 
theories have been discussed and evaluated by Walker^, and Chan et. al 1 ®. 

One example of a microphysi cally based constitutive law is an 
elasti c- viscoplastic theory based on two internal state variables as 
proposed by Bodner, et al. 11 . These authors, demonstrate the ability of the 
constitutive equations to represent the principal features of cyclic 
loading behavior including softening upon stress reversal, cyclic hardening 
or softening, cyclic saturation, cyclic relaxation, and cyclic creep. One 
limitation of the formulation though is that the computed stress-strain 
curves are independent of the strain amplitude and therefore too "flat" or 
"square" . 

Miller 12 h a3 reported research on the modelling of cyclic plasticity 
with "unified" constitutive equations. He also recognizes the shortcomings 
of many theories in predicting hysteresis loops, which are oversqaure in 
comparison to observed experimental behavior. Improvement is accomplished 
by making the kinematic work-hardening coefficient depend on the back 
stress and the sign of the nonelastic strain term. Theories that are 
similar in format to Miller's have been proposed by Krieg, Swearengen and 
Rhode 1 3 and by Hart 11 *. The models use two internal state variables to 
reflect current microstructure state and are based upon models for 
dislocation processes in pure metals. All these constitutive theories were 
formulated without the use of a yield criterion. Since these models do not 
contain a completely elastic regime, the function that describes the 
inelastic strain rate should be such that the inelastic strain rate is very 
small for low stress levels. Theories with a yield function and a full 
elastic regime have been developed for the case of isotropic hardening by 
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RobinsonlS^ an( j by L,ee and Zavrel^ for both isotropic and directional 
hardening. 

As previously noted, the quantities utilized in the small strain 
theory of viscoplasticity (stress, strain, stress rate, and strain rate) 
are defined only within the assumption of "small strains". Yet the precise 
definition of what constitutes "small strain" is always left unstated. 
Whether or not the stresses for a given case are "small" cannot be 
determined a priori by geometric considerations. In general, one cannot 
know in advance whether for a given loading of a material the "small 
strain" assumption (always left undefined) will hold or not. The question 
of whether the small-strain approximations are valid is always avoided in 
the "small strain" literature. Furthermore, as Hill^? points out, the 
really typical plastic problems involve changes in geometry that cannot be 
disregarded. In many cases, for example, it is sufficient to take into 
account finite plastic strains and small elastic strains or vice versa. 
From the theoretical viewpoint it is desirable in all cases to have a 
theory which intrinsically allows for both the elastic and plastic strains 
to be large. Such a theory of course, must reduce to the earlier mentioned 
special cases, as limiting cases. Furthermore, such theories provide a 
check for those which are obtained by generalizing small strain theories. 

The mathematical theories of deformation and flow of matter deal 
essentially with the gross properties of a medium. Heat and mechanical 
work are considered as additional causes for a change of the state of the 
medium. The resulting phenomena in any particular material are not 
unrelated. Therefore, a thermodynamical treatment of the foundation of the 
theory of flow and deformation is appropriate, and indeed the obvious 
approach. Two very different main approaches to a thermodynamic theory of 
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a continuum can be identified. These differ- from each other in the 


fundamental postulates upon which the theories are based. An essential 
controversy (a good survey of this controversy is given in Ref. 18) can be 
traced through the whole discussion of the thermodynamic aspects of 
continuum mechanics. None of these approaches is concerned with the atomic 
structure of the material. They, therefore, represent purely 
phenomenological approximations. Both theories are characterized by the 
same fundamental requirement that the results should be obtained without 
having recourse to statical or kinetic methods. 

Within each of these approaches there are two distinct methods of 
describing history and dissipative effects: the functional theory 1 9, in 

which all dependent variables are assumed to depend on the entire history 
of the independent variables, and the internal variable approach 20 , wherein 
history dependence is postulated to appear implicitly in a set of Internal 
variables. For experimental as well as analytical reasons 21 * 22 the use of 
internal variables in modeling inelastic solids is gaining widespread 
usage, in current research. The main differences among the various modern 
theories lie in the choice of these internal variables. 

The predictive value of an elastic” viscoplastic material model for 
non- isothermal , large deformation analyses depends therefore on three basic 
elements: 

a) the nonlinear kinematic description of the elastic-plastic deformation. 

b) thermodynamic considerations 

c) the choice of external and internal thermodynamic variables. 
a 3 well as on their interactions. 

The problem of viscoplastic deformations in shells has been treated at 
several levels of approximation and generality. 
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The simplest approaches* are based on the assumption of infinitesimal 
displacement gradients (which implies infinitesmal strains) and a material 
model of stationary creep, sometimes with an approximate inclusion of 
primary creep. 

A more general analysis utilizes shell kinematics for moderately large 
displacement gradients (at least some of them), infinitesimal strains, and 
material models of stationary or simple non- stationary creep . Ths type of 
assumption is capable of solving problems of creep buckling^, and it does 
reproduce the sometimes stiffening effect of the interaction between the 
normal forces and the normal deflection. Extension of these kind of 
formulation with a viscoplastic material model is presented in Refs. 25, 
26, it 27. The use of numerical methods^® makes possile the solution for 
many non- trivial types of structures. 

The problems of large strains, which arise in the analysis of large 
creep or thermal deformation of shells, have not been treated at all in a 
general manner. Recognizing that finite strain effects are present in 
these problems, reliable prediction demand tht such effects be included 
rationally and properly in the analysis. In addition to the necessary 
basic kinemtical and dynamical equations of the shell theory, such an 
analysis must incorporate a correctly invariant formulation of the material 
equations and requires an evaluation of the strain— rate tensors through the 
thickness of the shell. Such an analysis cannot be found in explicit form, 
at least in the readily accessible engineering literature. 

Several authors have developed mathematical description of the 
kinemation of the three dimensional deformation of elastic or 


* A comprehensive survey of these works is given in Ref. 23. 
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viscoelastoplastic materials^ » 30 # However, it is not clear how to best 
select to reference space and configuration for the stress tensor, bearing 
in mind the rheologies of realistic materials. Although an intrinsic 
relation, which satisfies material objectivity can be used3 1 »32 > the 
choice is not unique (see for example Refs. 30, 33 > 3 * 0 » 

1.2 Purpose of the Present Study 

The objectives of the present research are to develop a general 
mathematical model and solution methodologies for analyzing the structural 
response of thin, metallic shell-type structures under large transient, 
cyclic or static therraoraechanical loads. Among the system responses, which 
are associated with these loads and conditions, are thermal buckling, creep 
buckling and ratchetting. Thus, geometric and material type nonlinearities 
(of high order) can be anticipated and must be considered in the 
development of the mathematical model. Futhermore, this must also be 
accommodated in the solution procedures. The results obtained from this 
analysis are compared with the available experimental data, as well as with 
results obtained from "small strain" analyses in order to ascertain the 
range of validity of the "small strain" approximation. 

1.3 Synopsis of the Present Study 

Secton 2 contains the concepts that are necessary for the development 
of a general "finite strain" theory for thin bodies with path-dependent, 
time- dependent and temperature-dependent .material nonlinearities. 

A complete true abinito rate theory of kinematics and kinetics for 
continuum, without any restriction on the magnitude of the strains or the 
deformations is formulated. The time dependence and large strain behavior 
are incorporated through the introduction of the time rates of the metric 
and curvature in two coordinate systems; a fixed (spatial), and a convected 
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(material) one. The details of the reported developments are carried out 
by using tensor analysis. Special attention is directed to coordinate 
transformations as applied to continuum mechanics. Consideration of the 
kinematics of space (i.e., the intrinsic rates of change as they are 
observed by a geometer within a closed neighborhood of material particles) 
focused on the distinction between a fixed (spatial) coordinate system and 
a convected (material) coordinate system. Moreover, the rates of change of 
tensors are presented in both systems, taking the various tensor and 
kineraatical effects into account. The relations between the time 
derivative and the covariant derivatives (gradient) are developed, and 
these illustrate the possibilities of curved space or motion. 

The time derivatives are applied to the basic laws of continuum 
mechanics and thermodynamics to generate the equations of equilibrium rate 
and those for the compatiblity of the deformation rates. Finally, the 
principles of the rate of virtual power and the rate of conservation of 
energy are introduced and employed, and these should provide a basis the 
development of computational methods. 

The general form of the constitutive equations, employed in the 
analysis, is presented in Section 3. 

The metric tensor (time rate of change) in the convected material 
coordinate system is linearly decomposed into elastic and plastic parts. 
In this formulation, a yield function is assumed, which is dependent on the 
rate of change of stress, metric, temperature, and a set of internal 
variables. Moreover, a hypo-elastic law is chosen to describe the 
thermo- el as tic part of the deformation. 

A time and temperature dependent viscoplasticity model is formulated, 
in this convected material system, to account for finite strains and 
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rotations. The history and temperature dependence are incorporated through 
the introduction of internal variables. The choice of these variables, as 
well as their evolution, is motivated by phenomenological 
thermodynamic considerations. The nonisothermal elastic-viscoplasti c 
deformation process is described completely by "thermodynamic state" 
equations. 

Many pitfalls in the analyses of various investigations are indicated. 
A very important point that has been consistently neglected by many 
analystsj and computer programs is to indicate precisely in what form the 

I 

constitutive properties have to be input. Most investigators after an 
elaborate treatment of a general theory in tensor notation, leave undefined 
the constitutive equations to be measured in the laboratory. In Subsection 
4.1, the homogeneous uniaxial irrotational deformation of a continuum is 
treated, with at least two purposes in mind: (1) to give a clear physical 

understanding of the quantities involved in the analysis (which is not 
possible to obtain through the tensor index notation) and (2) since the 
most common material test is the uniaxial test, to identify precisely what 
are the quantities that one should measure in the laboratory (as well as 
how to express these data to conform with the constitutive equations used 
in the theoretical material model). 

In Sections 5 and 6 , the previous developments of Sections 2 and 3 are 
utilized to derive consistent kinematic, dynamic and constitutive 
equations, which are valid for finite strains and rotations of thin bodies. 
Some of these equations seem to be original (have not been found in the 
literature by the authors). 

Five different shell theories (approximations) in rate form, starting 
with the simple Kir chhof f-Love theory and finishing with a complete 
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unrestricted one, are considered in Section 5. Three different curved and 


straight beam problems are studied in Section 6. The results from the 
"finite strain" analysis are compared with the results from the 
"small strain" theory to ascertain the range of validity of "small strain" 
theory for the present type of problems. Moreover, the time and 
temperature dependence and effects of the new constitutive relations are 
compared with the results of the classical formulations of 
thermo-elasto-plasticity . 

The future research is summerized in Section 7. 
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2. GENERAL FORMULATION; THREE-DIMENSIONAL CONTINUUM 


2.1 Introduction 

The inherent difficulties associated with nonlinear continuum 
mechanics, along with theoretical considerations, lead to the formulation 
of incremental constitutive equations. Thus, the rate concepts arise and 
reformulation of the basic laws with the aid of the rate approach becomes 
essential. As usage of tensor operators is a common practice in the scope 
of continuum mechanics, their rates need to be formulated rigorously. 

It is obvious that the subtleties of tensor analysis and of the rate 
concepts, distinguishing between the spatial and the material descriptions 
of the continuum, require careful examination and rigorous formulations; 
otherwise, inaccuracies are likely to be encountered, especially in 
conjunction with the derivation of the rate of a gradient. 

The objectives of this section are then the systematic formulation of 
the rates of tensor operators, yielding integral and differential rate 
theorems. Resultant theorems are definitely meaningful. Moreover, the 
consistent considerations throughout the process of formulation are vital 
ingredients for gaining insight, in anticipation of further developments, 
when dealing with special structural elements. 

Treatises that have influenced this write-up are attributed to 
Truesdell et al.35“37 f Sedov38“ 1 * 0 , Green et al.^ 1 "^, Sokolnikoff 1 * 1, and 
McConnell ^5. 

The tasks of the reported developments are carried out with the aid 
of tensor analysis. Tensor definitions, notations, theorems, corollaries, 
etc. are rephrased in Subsection 2.2. This brief refresher on tensor 
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analysis follows Sokolnikoff* 4 * 4 and McConnell* 45 . Special attention is drawn 
to transformation of coordinates, bearing in mind the applications of 
continuum mechanics. 

Subsection 2.3 is dedicated to the kinematics of space, i.e., the 
intrinsic rates of change, as they are observed by a geometer within the 
closed neighborhood of material particles. The distinction between 
coordinate systems (essential cornerstone of tensor analysis) is 
specialized to the one between the fixed (spatial) system and the convected 
(material) one. Only the Lagrangian description is discussed because the 
Euler i an one does not seem adequate for solids. The geometric rates are 
examined in detail, as long as the geometer, whose explorations focus upon 
the material system, is capable of observing. This subject has been 
thoroughly elaborated upon by Durban and Baruch* 45 , whose discipline follows 
the original ones of Gibbs 14 7 or Bloch* 4 ®, and by Mendelssohn and Baruch* 49 
whose discipline in turn follows the work of Aris 50 . The latter two 
illuminate the facets of the various time derivatives (namely, rates of 
change) with the help of the conventional tensor notation. This approach 
will be followed, herein. 

Subsection 2.4 consists of the development of the rates of tensor 
components. Numerous publications (including textbooks), especially those 
oriented toward fluid mechanics, have dealt with the problem of rate of 
change of tensor components. Since the present work is concerned with 
solid mechanics, the various derivatives are going to be associated with 
material points and not with the spatial ones. Unlike the conventions of 
fluid mechanics (cf. Aris 50 ) and following the ideas summarized by Durban 
and Baruch* 46 , and by Mendelssohn and Baruch* 49 , rates of change must be 
distinguished by the location of the observer. 
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The intrinsic description is insufficient, as some basic laws refer 
to rates, which are observed within the fixed system. The partial time 
derivative and the total derivatives, although defined with respect to the 
fixed system, are formulated in the material system in a manner that 
manifests their tensor character. Besides, rigorous formulation of the 
time derivative of a covariant derivative (gradient) enables the 
classification of problems with regard to the curvature of space or motion. 
The time derivatives are applied to the basic laws of continuum mechanics, 
to obtain the equtions of equilibrium rate and the compatibility of 
deformation rates. Finally, the principle of the rate of virtual power is 
derived. 

2.2 Geometry and Statics of the Continuum 

Let the n-dimensional space be decribed by two systems of coordinates 

x 1 and u a (Fig. 2.1) 



Point P can be characterized either by its own coordinates x*(P) or 

by u a (P). The same holds true for the neighboring point Q. As for the arc 
PQ, it is characterized by the differences Ax* and Au°. When Q approaches 
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P t differences are replaced by differentials' and the following 
transf ormation stems 


. i 9x* . a 
dx - du 

9u° 


. ct 9u . i 
du - — 7— dx 

9x 


( 2 . 2 . 1 ) 


1 n 

where dx and du are the components of the differential of the position 
vector. In Eq. (2.2.1), and from now on, the summation convention is 
implemented, i.e., a repeated (dummy) index is summed over its entire 
permissible range (from 1 to n). Any set of n components, assuming the 

values a a when observed in the system (u^ ,u n ) and the values A 1 when 

observed in the other system, and which obeys a transformation law similar 

to Eq. (2.2.1), 


A 1 - — 

9u° 


9u° .i 
— T A . 
9x 


( 2 . 2 . 2 ) 


is defined as "the contravariant components of a vector". Distinction 
between the systems of coordinates is made by using capital (small) letter 
and Latin (Greek) indices for quantities as they are observed in the 

1 G 

x -system and u -system, respectively. 

If there exists a scalar function <t> which is point dependent, then, 
its gradient obeys another transformation law, namely 


9<t> 

9u° 


9d> 9x 
9X 1 9u° 


Any set of n components a^ snf A^, 


_9u 

9x 


i a o 


9x 

9u° 


(2.2.3) 

obeying the transformation law 
A i# (2.2.4) 
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is defined as "the set of covariant components of a vector". It should be 

i. o 

noted that the components of the position vector (x or u ) often do not 
constitute components of a vector, according to the definitions of Eqs. 
(2.2.2) and(2.2. J 4). It is obvious that "superscript" indices denote 
contravar iant components and "subscript" indices covariant ones. Moreover, 
it can be shown that summation conventions are applicable in a "diagonal" 
manner . 

The metric properties of the space are determined by the length of 
the elementary arc (as Q approaches P; see Fig. 2.1) in differential form 


ds^ « G^jdx^dx^ (2.2.5) 

where ds is the arc length. The quadratic form, Eq. (2.2.5), is specified 
by the elements (G^j) of a symmetric positive definite matrix. If the 
space can be described by a system (say y^) of rectangular Cartesian 
coordinates, the form of Eq. (2.2.5) reduces to 


ds^ - dy^dy* « Sjjdy^dy^ (2.2.6) 

where 6ij are the Kronecker deltas. However, points are observed in the 
u°“system too, hence 


o r* ft 

ds^ - g „du a du p and 
a{5 


3x* dx^ G 
‘ 3u° du 8 1J 


(2.2.7) 
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The sets G and g are defined as the covariant components of the metric 
i j 08 

tensor. It is symmetric tensor of rank two. Any set of n 2 components, 
obeying the transformation rules 

.lj 3X 1 oB 

° 77 du B 

»b- — *4 •?« <2 - 28) 

9u a 9x J 


* is! is! , 

« “ ax 1 9x J 06 

is defined as the set of contravariant (mixed, convariant) components of a 

tensor of rank two. A tensor is symmetric if it is not affected by an 

interchange within a pair of indices (i.e.. Ay - Aji, A*J - aJ 1 , etc.). 

The dots in the mixed components help to distinguish between a* 8 , a^.etc. 

The conravariant components of the metric tensor are readily 
available by matrix inversion, namely 


G 1J G - 6 - , 
Jr r 


o$ -a 

8 8 fnr " 


( 2 . 2 . 9 ) 


It is obvious that the mixed components of the metric tensor are the 
Kronecker deltas. The metric tensor facilitates the following operations: 
raising indices 


a* ■ 


lowering Indices 


ft co * B pB 

a a * g a6 a 5 a a6 “ g a ? g Bc a ! V “ g ap 3 


( 2 . 2 . 10 ) 

( 2 . 2 . 11 ) 
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inner product 


•\ ' S°V 6 ' - •.*' (2 - 2 -' 2) 

No transformation rule is violated by Eqs. (2.2.10) and (2.2.11). 

From now on, the mixed components of a tensor are going to be witten 
without dots, as no raising or lowering is intended. The idea is to keep 
track of the omitted dots and remember that their loction is identical in 
any system of coordinates. Then, the set of nP^ components, obeying the 
transformation rule 


V--J, 


ax 


L i 


B, 


9u 


9u 6 q « r .. o p 

Jq V" *q 


(2.2.13) 


9u ax ax 

is defined as the set of the covariant (rank p) and covariant (rank q) com- 
ponents of a tensor. Raising or lowering an index affects the nature of 
the transformation, Eq. (2.2.13)* 

It often occurs that although a set of indexed quantities can be 
obtained, the verification of the tensor character of such a set by 
inspection (whether the set follows the transformation law of Eq. (2.2.13)) 
may be quite cumbersome. The inconvenience of the verf if i cation of the 
transformation laws can be circumvented with the aid of the quotient 
rule, as follows: let A(i ^...ip, j r *..J q ) denote a set of nP+Q quantit- 

ies, the indices of which are ^ through i p and J 1 through J q , and let 

denote the components (contravariant of rank s and covarariant 

i 1 i r 

of rank r, (r < p, s < q) of an arbitrary tensor, if 


j 1 • • • J , 


(i r .., i p ; j 1 ..*j q ) Bj A 


i. . . 

, 1 pr 

'j ... 


( 2 . 2. 1 *0 
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where the C's in Eq. (2.2.14) are the components of a tensor, and the A’s 

i 1 i 

are also components of a tensor, namely A. This rule is very 

J 1 J 1 

useful in the course of tensor algebra and tensor calculus as relations 
such as Eqs. (2.2.14) are often by-products of various derivations. Thus, 
tensor characteristics can be easily verified (or identified). 

It is often required to calculate n-tiple integrals over a specified 
domain in the n-dimensional space. The n-tiple "volume" element is 
calculated as follows: let the n-tiple parallelopiped be constructed from 

a vertex at point P, the edges along the coordinate lines (dx 1 for the i-th 
edge) and the opposite vertex at point Q, the coordinates of which are 


obtained by adding dx 1 to the coodinates of P: then, the volume dV is 

dV - Vg* dx\.. dx n (2.2.15) 

where G is the determinant of the square matrix, the elements of which are 

the covariant components of the metric tensor. 

Christoffel symbols of the second kind are defined by the formula 


, « 3G. 

1 G is f — J2 ♦ 

2 1 k 1 

* 3x K 3 x j 


3G. 3G,. 

ks jk 


) 


J k 9x k 

These symbols are not tensor components. However, 

transformation rule 


( 2 . 2 . 16 ) 

Christoffel derived a 


3 2 x r 

3u a u® 


,p M_ 
06 3u P 


- r 


3x J 3x 
Jk 3u a 3u e 


(2.2.17) 


where the Ps denote the Christoffel symbols, derived from the g's in the 


u a -system according to Eq. (2.2.16). It is obvious that Christoffel symbols 
do not transform like tensor components unless the coordinate transforma- 
tion (between the u's and the x's) is affine. 
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There are several tensor operations, in order to produce tensors, 
e.g,: addition (tensors of same rank), outer products (without repeated 

indices), contractions (as implied in Eq. 2.2.14)), etc.. All these 
operations are algebraic in nature. However, tensor fields (i.e., the 
components depend on the coordinates) can be operated upon 
bydifferentiation in order to produce tensors of higher ranks. It can be 
shown that the covariant derivative of a tensor, defined as follows 
i,...i * i, . . . i, P I. i....r...i_ 


1 * * p _ _J_ A 1 * * P + j , * A T” 
i 1” ,J q,k 9x k J T'* J q A-1 P k V” J < 


? r r , 1 1 * * * 1 p 

' VV"r-J, <2 - 2 ' ,8) 

does actually transform according to Eq. (2.2.13) as a set of components 
contravari ant of rank p and covariant of rank q + 1. The comma in Eq. 
(2.2.18) denotes the covariant derivatives not to be confused with partial 
derivatives with respect to the spatial coordinates. The combinmation of 
dots and of the index r in Eq. (2.2.18) means that the repeated (dummy) 
index replaces the Ith subscript (or superscript), for the sake of contrac- 
tion. It should be noted that the covariant derivatives of a tensor 
really show how It varies while the observer is moving along a coordinate 
line. Hence, the operators of vector analysis (gradient, divergence, curl, 
etc.) have to be phrased in terms of covariant derivatives. 

As a consequence of the definition, given by Eq. (2.2.18), it can be 
verified that covariant differentiation of sums and products should be 
accomplished in the same manner as partial differentiation does. Moreover, 
Ricci’s lemma 


g o$,v - 8 ’ v “ 0 


(2.2.19) 
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suggests that the components of the metric tensor behave like constants as 
long as covariant differentiation is concerned.. Thus, the operations of 
raising and lowering indices and (on the other hand) of covariant 
differentation are commutative. 

Naturally, covariant derivatives of higher order are obtained 
systematically, i.e., 


i .... 1 
\ J P 

J T ,,h q,kl 


i, ...i 

U J j' ] 
J T’’ J, q,k ,1 


( 2 . 2 . 20 ) 


The order of differentation may be significant, as 


P 

A . . , - A . , . - R* .. A 

i,jk i#kj ijk p 


( 2 . 2.21 ) 


where the R*s are the component of the tensor of Rieraann-Christoffel namely 


r p - _L r p + r 3 r p - r s r p 

ill. Ml. 1 J 4 4 1 


•ijk ax J ' ik 3x k 1J *ik *sj ij 3k 


( 2 . 2 . 22 ) 


This tensor manifests the curvatures of space. If and only if, the 
Riemann-Christoffel tensor components vanish identically, then the space is 
Euclidean (i.e., can be represented by a system of Cartesian coordinates). 


The transformation tensor x is defined by 

a 

i A 3X 1 


( 2 . 2 . 23 ) 


It can be proven that the quantities x q transform like the contravariant 
components of a vector in the Latin indices (say, if the x-coordinates 

transform to another system of x-coordinates) and like the covariant 


components of a vector in the Greek idices (if the u-coordinates transform 
to another material system of u-coordinates). Therefore, the symbol x^ is 
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written (following Sokolnikoff Ref. [44]) with a Latin superscript and a 
Greek subscript. A similar reasoning yields the definition of the inverse 
transformation tensor 


a ^ 3u a 

. i 
9x 


(2.2.24) 


A treatment of double tensor, such as x^, requires a definition of 
total covariant derivative or the tensor derivative. Let a double tensor 
be given by 

A®***? (x 1 ^ 0 ) (2.2.25) 

T»»«0 k«*«i 

Its total covariant derivative is defined by 


■ Qt • . . S i . • . J ^ » o . . . 0 i . • « J ^ .a . . . 0 i • . . J # ^ 

•y...6 k. . .1; r ’ y...6 k...l,r >...« k l,m x f 


where, A 


ct . . . 8 i . . . j 
y . . . 6 k . . . 1 , f 


9A 


O...0 1...J 

Y, ‘- 6 — — — ♦ y a A?"*? , i ***'J + ...(2.2.27) 

f pf >.,.o k...l 


3u 

_ P + 

Sr p...6 k. • ,1 

Thus, it is a partial covariant derivative with respect to u°, where x* are 
constant. For tensors which are not double, the definition given by Eq. 
(2.2.26) reduces to the regular one. 

i i Q 

The transformation tensor x fulfills (x depends on u only), 

ct 


3x 

9u 2 


i 


1 A i 
X ;a “ *a 


(2.2.28) 


The tensor formulation of the basic equations of continuum mechanics 
has been presented in numerous books, e.g.. Fungal, Sokolnikoff^* Green 
and Zerna^l f etc. Let the three-dimensional Euclidean space (i.e., our 
empirical engineering space) be described by the system of coordinates x^- . 
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The solid continuum occupies a definite volume in space, the volume element 
is denoted by dV. The continuum is bounded by a surface, the area element 
is denoted by dA, v* denote the contravariant components of the unit normal 
to the surface, i.e., a unit vector that is perpendicular to the surface, 
namely 

G ljV i v J - 1 (2.2.29) 

An area element consists of specific material particles. Let v stand 
for their density. This concept is analogous to that of mass density, 
however, it is intended to specify a measure for the material entities and 
not the mass of inertial considerations. Since no dynamics are involved, 
let v be called the "area mass density" and let dm be the "mass of an area 
differential element, i.e., 

dm - v dA (2.2.30) 

If dP 1 denote the components of the force vector, acting on the area 
element dA, then the surface tractions are defined as the mass density of 
that force vector 

T 1 - • T 1 =* -£§7 ' YT 1 * (2.2.31) 

dm > TdA / dA 

The last equality of Eqs. (2.2.31) defines the "conventional" surface 
tractions. The present definition, Eqs. (2.2.31), is based solely on 
grounds of convenience in view of the foreseen kinamatic developments. 

The contravariant components of the stress tensor, as defined by 
Cauchy (cf. Fung 51 , Sokolnikoff ^ etc.) are denoted by a^-J, and they are 
associated with the tractions and with the unit normals as follows 

rr J « (2.2.32) 


2h 



Thus, it is obvious that the stress components are associated with a 
specific volume, whereas the components of the second P iola-Kirchhof f 
stress tensor (cf. Hibbit et al.^2) 

a iJ - 4 iJ (2.2.33) 
p o 

are associated with mass (material entities) as stated in conjunction with 
the traction, where p and p 0 stand for the current mass density and for a 
reference one, respectively. 

As long as the "classical" terminology of tensor analysis is 
concerned (cf. Sokolnikoff , the components of the stress tensor (aiJ) 
are those of absolute tensor, while those of the Piola—Kirchhoff stress 
tensor (S*-)) constitute a relative tensor. 

Let fJ denote the components of the body forces per unit mass. The 
equations of equilibrium are then 


o 1 ^ + pf J - 0 (2.2.34) 

accompanied by the boundary conditions, Gqs. (2.2.32) at the tracted 
boundaries . 

These equations are based on equilibrium considerations exclusively. 
The kinematic constraints are systematically incorporated by the principle 
of virtual power (several publications name it "The Principle of Virtual 
Velocity"). Let 5vj denote the components of the vector of virtual 
velocity, i.e., a vector field of values of velocity at which the material 
particles are capable of moving, obeying all kinematical restrictions 
(continuity, prescribed velocities where the virtual ones vanish, etc.). 
Then, the integral theorem 
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(2.2.35) 


j o ij 5Vj i dV - | pf J 6v^dV - j 7T J 5VjdA - 0 

V V A 

is equilvalent to the equations of equilibrium and to the complete set of 

boundary conditions (Note that 6vj = 0 implies that vj is prescribed). 

2.3, Kinematics of the Continuum 

Let a continuum be described by the two coordinate systems of Fig. 
2.1 The x^—system stays at rest and will be labeled as "the fixed system", 


while the u°-system is associated with the material points and will be 

labeled as "the material system". As the continuum is moving and 

deforming, coordinates x* of a material point (say P) are changing, i.e., 

x!(P) are time-dependent. Yet, a material point preserves its identity and 

hence its material coordinates (u ) are not changing in time (t). The 

o 

coordinate lines are assumed to "deform" with the continuum in order to 

enable the material points to keep their coordinates (in the u°-system) 
unchanged . 

The components of the velocity vector in the fixed system are defined 
by: 



While a geometer at P is observing, he cannot recognize any change in the 
coordinates of any point, however, the main purpose of a geometer’s 
interest, namely distances, are obviously changing. As a matter of fact, 
the differential form of Eq. (2.2.7) should be examined. To do so, 
following Ref. 35, we must define material derivatives. 


Let A 


i . . . . j ot . • . . 0 
k. . . . 1 Y....5 


A i....j a....B ^i^ u a^ be aome double tensor 

k. . . .1 Y. . . . o 
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DA* * * 

depending on both x 1 and u q . Then its material derivative is given 

by 


DA 

Dt 


i 

k 


j a 
i y 


B 

6 



• j a 
.1 y 


B 

6 


* ... 


S A 1 

r, A 
kq s 


. j a. . . 0 

.1 y...6 



const 

const 



+ 



j.i ft s. . . j a. . . 6 
sq k...l Tf,..{ 


( 2 . 3 . 2 ) 


+ A*** . v 1 

• • • »q v 


u a - con^ 
x 1 - const 

The general definition of Eq. (2.3.2) has two particular cases: 

A) When one uses only the fixed spatial description for the description of 
tensor A, i.e., A = A(x*,t) (as is common in fluid mechanics), then its 
material derivative, Eq. (2.3-2), becomes 


DA 


i 

k 



Dt 


J a- 

1 y. 


B 

S 


3A* ” 

♦ t I 

3t 


. ♦ A***' 

x - const • • • » Q 


(2.3.3) 


The second term in the right hand side of Eq. (2.3-3) is called the 

translation term of A”* . 

• * • 

B) When one uses only a material description for A, i.e., A - A(u a ,t) (as 
is common in continuum mechanics) the material derivative becomes 


. DA 
Dt 


i.. 
k. . 


. j a. 

.1 y. 


.0 

.6 



const 


i .s. . .J o. . . 8 
sq A k...l y...S 


- r 


* A 1 
kq s 


.J a. ..6 _ i q 

.1 y...6 ***' v 


(2.3.4) 
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In this case, ~ - - is also called an intrinsic derivative of A. It is 

U w 


worthwhile to notice that, for single tensors in the material description, 
the definition in Eq. (2.3.4) reduces to a simple partial derivative: 


DA 

Dt 


a 

y 


. 8 
.6 



u° - const. 


(2.3.5) 


For the velocity vector, v*, one can show that the definition given by Eq. 
(2.3*1) can be replaced in the following manner 


i dx* 9x* Dx* 

v “ dT ■ 3T “ dT (2 - 3 - 6) 

From here on, all tensors will be given in the material description only, 
unless otherwise specified. So all material derivatives of the different 
tensors will be given by Eq., (2.3*4) (for the most general case of double 
tensors ) . 

Differentiation of Eq. (2.2.7) yields 


fe (os 2 ) - 203 (d.) - 8 oB ) duW (2.3.7) 

By dividing both sides of Eq. (2.3*7) by 2ds 2 , one obtains 

§j(lo 8 d3).I(kg aS )3f £ 6 <2.3.8. 

. a 

Obviously, is a unit vector, and Eq. (2.3*8) expresses the logarithmic 
deformation rate. The clue for the intrinsic rates of change may be 
unraveled, then, by the derivation of the rates of change of the metric 
tensor. Hence, the second of Eqs. (2.2.7) should be differentiated to 
yield 
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(2.3.9) 


3t ^ g a0^ “ Dt ^ X a X 0 G ij^ 

As attention is drawn toward the differentiation of Eq. (2.3.9), it 

should be noticed that the components of the transformation tensor are 

associated with a material point (they are changing while the continuum is 

moving). Nevertheless, the components of the metric tensor are associated 

with a spatial point and their sole cause of variation is the variation of 

the x- coordinates. Hence, they are inherently independent of time. Then, 

time differentiation yields 


. Dx 1 , . Dx J 

3 _ r f - J + y 1 - 1 

9t g oB °ij 1 Dt 0 a Dt J 


( 2 . 3 . 10 ) 


The time derivatives of x^ follows the general form of Eq. (2.3.^) 


Dx 


Dt 


3x 

c 

at 


a 


3u 


. r 3x v ,.1 s k 

+ r - x - xv " — ( — ) + r sk V 


.i i s k a 
, x , x 
sk sk a 


dt ' au a 


i 


( i2L ) ♦ r i x s v k 

1 3t ; sk a 


(2.3.11) 


The last equality in Eq. (2.3.11) is Justified by the fact that material 

coordinates (u°) are time - independent. Use of Eq. (2.3*6) into Eq. 

(2.3.11) and of the definition of covariant derivatives, Eq. (2.2.26) 
yields 


Dx 

o 

Dt 


3,1, i k s i i 

Iv)+r.vx - V -V 

. a 1 sk a ;a ,a 

3u 


( 2 . 3 . 12 ) 


Since v^ denote components of a vector in the fixed system, they can be 
transformed to the material system by Eq. (2.2.2) 

1 1 ° -P • p -- r (2.3.13) 


i o 

X V 
a 


P p r 
v K » u v 
r 
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Naturally, the velocities in the material system do not measure time 
derivatives of coordinates; their sole significance stems from the 
transformation rule, Eq. (2.3.13), as the geometer cannot measure them 
directly (this argument will be dwelt upon later). Hence, thusly operating 
on Eq. (2.3.12) yields 


i .ip. i p^ip 

v - (x v K ) - x v K + x v r 


,o p,o 


p ,a 


Special treatment is required for the evaluation of x 


P, a 


(2.3. 1*0 
As already 


stated, x* are tensor components in both coordinates, p & a systems, then by 
Eq. (2.2.26) 

3x k 


t 1 x 1 - > p x 1 + r* x J 

a,8 9u 6 a 08 P Jk a ^6 


3 2 x l 


3 u° 3 u B 


- v p „x* + r* x J x k 


08 P jk a 8 

However, substitution of Eq. (2.2.17) into Eq. (2.3.15) yields 


(2.3.15) 



and finally 


( 2 . 3 . 16 ) 


D r i 1 1 ip 

rrr- X J - V - XV 

Dt a ,o p ,a 


(2.3.17) 


Now, Eq. (2.3.10) may be expressed in terms of velocities (deformation 
rates) 


g * G.. (x i v p x^ + x 1 x^ v p ] 

Dt ®a 8 ij p ,a 8 a p ,8 


o p 

B V + St V * V + V 

B p8 ,a B ap ,8 a, 8 8, a 


(2.3.18) 
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The components of the deformation rate are then, defined as 


oB 


— f v + v ) 

2 1 o,8 p,<* 


Bo 


(2.3.19) 


and thus, 


_ o . p 

Dr, j v . du du 

Dt (l ° g d3) ‘ d aB dl" di" 


( 2 . 3 . 20 ) 


Dt ^ g aB ^ " 2d aB 


( 2 . 3 . 21 ) 


The intrinsic geometer is, by no means, capable of measuring 
velocities directly, Eq. (2.3.13). yet he can measure metric rates, Eq. 
(2.3.18), and obtain the velocities by spatial integration of Eq. (2.3.18). 
Two problems arise in conjunction with this: (a) a system of equations 

like Eqs. (2.3.19) (i.e., given d and required v is not necessarily com- 

Op o 

patible, and (b) even if compatibility is ascertained, the solution is not 
unique. Treatment of the first problem will be postponed, while the second 
is alleviated by the definition of the spin tensor 


then 


”<,6 ■ 5 (, „,6 • T *.«1 • ' “eo 


V o,B " d aB * W a8 d Ba W 0a 


(2.3.22) 


(2.3.23) 


a 8 


The derivatives of the contravariant components g are readily 
obtained by the differentiation of Eq. (2.2.9), namely 

JL (g aS g ) 4- (6 a ) - 0 (2.3.24) 

3t ^ g g 0V at y 

hence 

(2.3.25) 


3 oB _ 0 oB . 

g BT at g 2g d B^ 
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( 2 . 3 . 26 ) 


oy 

Multiplication by g and contraction yield 

2 *° 6 s P \v 

As soon as the deformation rate is established as the time derivative 
of the metric tensor, the intrinsic characteristics of the continuum (i.e., 
Christoff el symbols, the Riemann-Chr istoffel tensor, etc.), being metric 
properties of space, are readily differentiated (with respect to time). 
This depends on the deformation rate (d^g) an< * hot on the spin (w^). As 
an illustration, the determinant of the metric tensor (g) is a measure for 
a volume element (a detailed proof appears in Ref. M9 and the result for 
the rate of dilitation is 

It c ■ G y !p (2 - 3 - 27) 

It is evident that the spin components do not affect the expressions 
in Bq. (2.3.27) (dilitation is not affected by the spin). 

As Christoffel symbols, Eq. (2.2.16), play an essential role in 
tensor calculus, the need for an evaluation of their time derivatives 
cannot be underestimated. Hence 


J_ yO . i f i_ ) r !f£i 4 - IffiS. ) 

3t pa 2 1 3t * ' 1 ^..c *„p ,..B 


3u 


3u p 


3u 


3g„ 0 


36. 


f __£i + - _££ ) 

O & ^4- ^ ~ #1 Q ' 


3t ' 3u° 


( 2 . 3 . 28 ) 


3u 


3u k 


Systematic treatment of the first term in Eq. (2.3.28) yields 


1 ( i K ° e ) . g - - 2^ W g° e d 

2 3t ® ‘ d po®Bp pd® Bp 


( 2 . 3 * 29 ) 
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Now since the u- coordinates are time- independent , then 


1 i_ r . i_ ( 1 ) . i_ d 

2 9t 3u° 3u" 2 at Su" (>8 


- d „ + y v d- ♦ y v B d (cf. Eq. 2.2.18) 

p8 ,a p a v6 Ba VP 

Interchanging the triads of indices and rearranging yield 


(2.3.30) 


i 08 9 r 9g pe . 9g 6o 

28 -k'Ts'ZF 


9 g 


Ba 


) - g 06 [d 


9u 


pB.a 


+ d„ - d . 
Ba,p pa, 8 

♦ 2 V v d J 
pa vB 


(2.. 3. 31 ) 


Substitution of Eqs. (2.3.29) and (2.3.31 ) into Eq. (2.3.28) yields 


9a aB r ^ . 'j 

3t Y pa “ 8 d p8,a + d Ba,8 ~ d pa,8' 


(2.3.32) 


In cases where rates like those in Eq. (2.3.32) are desired in terms of 
velocities, elementary manipulations lead to 


9 1 


» 1 r v a 4 v° 4 1 V [r e ° + r 5 a ) 

3t pa 2 1 ,po ,ap 2 8 *p*v *a*p' 


(2.3.33) 

It should be noticed that although Christoffel symbols are not tensor 
componets, their time derivatives are. Specialization of Eqs. (2.3*32) and 
(2.3.33) for Euclidean spaces are self-evident. 

The treatment of curvature rates concludes the discussion of the 
kinematics of the continuum. Regardless of the type of space considered, 
Eq. (2.2.22) should be differentiated. Then 


3 r a i _9_ r _3 yO 

3t r • pax^ 9t -a px 


° 4 y X y° - — y a - v A y a ) (2.3.3 1 *) 

3u v - pT >0 3 u t pa p0 * T 
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Changing the order of differentiation results in spatial derivatives or a 
tensor , namely 


JL r JL y a ) . i_ r JL y a ) 

9 p° PT 9 p° 9t pT 


(2.3.35) 


Benefitting from the rules of covariant differentiation 


9Y a 3V> 

9 r _ 9 _ t t PJ i _ y a pr 

. a 9t p t 9t ,o Ap 3t 

<7U 


+ y 


i 


9y 


It 


♦ y 


pp 9t 

The first pair of terras in Eq 



pt 3t ’ 

(2.3.34) reduces to 


(2.3.36) 


4- ( i~ ♦ y a y® 

9t g u P PT pT Ap 


9Y 

)- (iF) 


fO 


+ 7 


pa 


37® . 97® 

At + > 

3t a t 3t 


♦ Y A (2.3.37) 

PT 9t 

The second pair of terms is obtained by changing indices. Substitution 
into Eq. (2. 3. 3*0 and rearrangement yield 


9_ 

3t 



) - ( 







(d - d . + d. ♦ d . - d. ■ ) (2.3.38) 

PT.TP pA,PT iT,pP PT,Ap Ap.pT 


As a summary of the kinematic considerations, Eq. (2.3.28) is 
bivalent. On one hand it relates velocities to the logarithmic rate of 
deformation. This relation is very useful in applications where velocities 
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are the principal dependent variables. On the other hand, if the 
deformation rates are given, Eq. (2.3.23) changes its role from formula to 
equation, and it must be solved for the velocities. In the latter case, 
Eq. (2.3.38) is an essential asset, because it expresses the compatibility 

conditions, i.e., a symmetric tensor d^ can be introduced into Eq. 

(2.3.18) to yield the resultant velocities v^, if they (d^) satisty Eq. 

( 2 . 3 . 38 ) . The commonplace Euclidean space supplies an immediate 
manifestation. The the t’s and their time derivatives identically vanish, 
and the familiar equations of compatibility from the mathematical theory of 
elasticity (cf. Sokolnikoff Ref. [44]) appear. Kinematics and deformations 
of surfaces present another example, then an intrinsic relation like Eq. 

(2.3.38) must be equated to the known time derivatives of the 
Rieamann-Chr istof f el tensor (say, from the equation of Gauss, cf. 
Sokolnikoff^ 1 *) . 

2.4 Time Derivatives of Tensor Components 

Let a continuum in space be moving (translating, rotating, deforming) 
together with tensor fields associated with its material points. The 
meaning is that various tensor fields (e.g., forces, stresses, heat fluxes, 
etc.) are observed and measured in conjunction with material points. This 
is the Lagrangean viewpoint. However, the position of the observer is 
very significant. If the observer resides somewhere in the material system 
it relates the tensor components to a deforming system of coordinates. 
Thus, the rates by which the components change may appear insufficient for 
the complete understanding of the observed phenomena. On the other hand, a 
fixed observer (naturally, in the fixed system), while observing the 
various fields in conjunction with a specific material point (he is capable 
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of tracing a moving point) associates the components with time- independent 
metric properties. It should be stressed that the latter observer does not 
pretend to be Euler ian, it is Lagrangean as the former because its scope of 
interest (material points rather than spatial points) matters and not his 
residence . 

g 

Hence, let the symbol ^ denote the operation of material deriva- 
tive, by which the rate of change is measured within the material system. 
Consequently, the tensor components are observed in order to evaluate their 
rates of change; however, convective terms (cf. Aris^O) originating from 
the curvature of the coordinate lines, are also taken into account. On the 

other hand, the symbol — ~ is used to denote the operation of total deriva- 

QU 

tive, by which, the rate of change is measured from a fixed standpoint. 
The components, the rates of change of which are evaluated by the operation 
of the total derivative, are refereed to unchanging coordinate lines. As 
rates within such a frame of reference seem to be most adequate for several 
applications of mathematical physics (cf. Sokolnikoff 1 * 1 *, the rules of total 
differentiation have to be elaborated upon). 

Let, then, the transformation rule of Eq. (2.2.13) be rephrased with 
the aid of the definitions of Eqs. (2.2.23) and (2.2.24), namely 



(2.4.1) 


The total derivative is defined as the rate of change of the components 
i . .... 1 

A P (including the convective effect and those of curved coordinate 

Jl . . . .3^ 


lines); however, formal implementation of this definition on Eq. (2.4.1) 


leads to: 
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(2.4.2) 
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for the fact is that the partial derivatives of all the bracketed terms in 
Eq. (2.4.2) have already been designed to account for the convective 
effects. A doubt arises whether the expressions in Eq. (2.4.2) do actually 
yield a tensor out of a tensor. As a matter of fact, the sole "trouble" is 


embodied in the partial derivative of the tensor components. Let, then u*' 
denote another set of material coordinates, and let barred symbols denote 
quantities in the u a -system, then, following Eq. (2.2.13) 
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and the sums in Eqs. (2.4.4) vanish because the material coordinates are 
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inherently time- independent . Equation (2.^.4) expresses, then, a tensor 
transf ormation rule. 

Returning to the definition of Eq. (2.H. 2), formal operations imply 
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The derivatives in the first term of Eq. (2.4.5) are readily obtained 
from Eq. (2.3.17), those of the second sum may be obtained by 
differentiation of the identity 
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Multiplication by u B , contraction and rearranging yield 
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The terms in Eq. (2.4.5) should, then, be elaborated on 


(2.4.8) 



a 


} a 6 


1 ”* 
1 ’ " 





(2.4.9) 


38 



> , 8 », V"°p 

3t "V 6 r ..B q 


6 , 


1 P 

-V U . a. 

,P .J A 6 1 


V 0 p B i P 0 r- ,a p 

. B U ,J. V ,e®P 1 ...P...B, (2.4.10) 

q a h l q 


The latter equalities in Eqs. (2.4.9) and (2.4.10) stem from interchanges 
between repeated indices. Substitution of Eqs. (2.4.9) and (2.4.10) into 
Eq. (2.4.5) yields 
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This is a typical manifestation of the tensor character of the total 
derivatives. It should then be rephrased as the rule for total 
differentiation in the material system, namely 
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It is obvious that the total derivatives depend on the velocity gradients. 
Looking forward towards the influences of the various types of motion, Eq. 
(2.3.23) is utilized to obtain 
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where the symbol — denotes the operation of the Jaumann time derivative 
dt 

(cf. Fungal), i.e., total derivative with the rotational effects discarded, 
namely. 
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The total time derivative and the Jaumann one have two interesting char- 
acteristics: sums and products of tensors are differentiated with the aid 

of the same rules of the differential calculus, and the derivatvies of the 
metric tensor vanish identically. The former statement establishes the 
formal readiness of the time differentiation of tensor equations, while the 
latter enables lowering (or raising) indices before (or after) the time 
differentiation without affecting the form of the equations. 

Last but not least, the total dif f erentation of the covariant 
derivatives of tensors should be studied. The exact relation between both 
differentiations (temporal vs. spatial) is the clue for a reliable rate 
theory (and applications). Obviously, the covariant derivative of a tensor 
is a tensor, the covariant rank of which is higher (by one) from that of 
the original one, hence from Eq. (2.4.12). 
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The first term in Eq. (2.4.15) is not familiar, then by Eq. (2.2.27) 
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Now, the first term in Eq. (2.4.16) is obtained by the rule of total 
differentiation of Eq. (2.4.12), namely 
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Collecting terms of Eqs. (2.4.15, 2.4.16, 2.4.17) yields 
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Substitution of Eq. (2.3.33) and utilization of the characteristics of the 
Riemann-Christoffel tensor (Sokolnikoff 1 * 1 * presents them as exercises) yield 
the final result 
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In summary, the operations of time differentiation that of covariant one 


are permissible only in case of rigid body translations (v w - 0), and any 
other motion implies further terms. In a Euclidean space, the effect of 
higher tensor rank is added. Curved spaces imply more complicated terms 
originating from compatibility considerations. 

2.5 The Principle of the Rate of the Virtual Power 

The principle of virtual power (or of virtual velocity), as expressed 
in Eq. (2.2.35), is equivalent to a set of differential equations 
(equilibrium) and associated boundary conditions .(should be appropriately 
selected between kinematical and natural). However, practical reasoning 
concerning solution procedures gives rise to basic doubts, as follows. Let 
a system of specific body forces (fJ) and specific surface tractions (tJ) 
act on a solid continuum. It is, natually, deformed and stressed. The 
stresses satisfy Eqs (2.2.35) or (2.3.34) but their relations with the 
deformation are not self-evident. There are constitutive theories 
asserting that considerations of convenience (hyperelasticity) or necessity 
(hypoelasticity) introduce incremental "stress-strain" relationships, i.e., 
the stress rate is assumed to depend on the deformation rate. In such 
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cases, equations like Eq. (2.2.35) need be manipulated further in order to 
obtain terms with stress rates. 

The objective of this section is, then, within the scope of obtaining 
the total time derivative of the principle of virtual power. The resultant 
virtual theorem may be applicable for the derivation of rate equations 
(equilibrium rate and incremental boundary conditions) and for the direct 
development of solution methods (e.g., incremental finite elements). 
Throughout this section, the space is assumed to be Euclidean and all 
tensor components, following the notations at the end of Subsection 2.2 
(Eqs. (2.2.17) and so forth), are measured only in the material system of 
coordinates. Hence, the tensor notation reverts to Latin indices, keeping 
the material system in mind. 

Eq. (2.2.35) should, first, be rephrased 
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Comparison between Eqs. (2.5.1) and (2.2.35) may explain the symbol dm, it 
stands for the mass of a (material) volume element or surface element 
respectively. The volume and the area are varying during the deformation, 
the mass remains constant. Total diff erentiaton of Eq. (2.5.1) yields 
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The first integral in Eq. (2.5.2) needs detailed treatment, 
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The last term in Eq. (2.5.4) symbolizes the volume changes, however from 
Eqs. (2.3.21) and (2.3.27) one can obtain 
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The first term in Eq. (2.5.4) consists of the desired stress rates, and 
therefore the incremental constitutive equations should be substituted 
there. The remaining integrals in Eq. (2.5.4) should, then, be elaborated 
upon, since the time derivative of the (virtual) velocity exists. This is 
a typical need of the rate of a gradient, by Eq. (2.4.19) 
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Collecting terms and rearranging yield 
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At any instant, Eq. (2.5.7) is satisfied. The virtual velocity and its 

time derivative are, then, independent. Moreover, the last three terms of 

Eq. (2.5.7) (those depending on the "virtual acceleration") are equivalent 

to Eq. (2.2.35). Hence, the principle of the rate of virtual power may be 
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obtained in its concise form. 



do 

dt 


iJ 


- 1J V - ° kJ \*> 


6 v. . dV 
J I 1 


r 

J 

V 


df J 

p i- 6v jd V - 



6VjdA 


0 


(2.5.8) 


For further classifications, the total derivative of the stress components 
will be represented by Jaumann derivative, Eq. (2. 4. 1 1 *), namely 
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Let the following integrals be defined 
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Then, substitutions in Eq. (2.5.8) yield 
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The integral I e , Eq. (2.5.10), is expected to be formally similar to 
the "conventional" integrals of linear elasticity. The integral I r , Eq. 
(2.5.12), is the sole outcome of rigid rotations, and the integral Iq, Eq. 
(2.5.11), is required for completeness. 

It is obvious that further manipulations on Eq. (2.5.8), especially 
the use of Gauss* theorem, lead to the field equations and associated 
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boundary conditions. However, the integral form is quite satisfactory for 
the utilization of the finite element method. Further inspection of Eq. 
(2.5.8) should reveal the resemblance and the differences to the integral 
incremental principles presented by Hibbitt et al^ 2 and numerous successors. 
Their scope is quite different (the P iola-Kirchhof f stress tensor) and 
their formulation contains several simplifications. 

2.6 Concluding Remarks 

The formulations of this section focus upon tensor rate operators 
associated with continuum mechanics. Tensor definitions and corollaries 
are compiled. The deformation rates and related topics regarding the 
kinematics of space are presented. Compatibility of deformation is 
elaborated on as an outcome of the curvature of space. 

Beside the intrinsic rates of change, total (time) derivatives and 
partial derivatives of an arbitrary tensor are presented. The material 
system of coordinates is useful for applications, the spatial is favorable 
for reference, their correlations are highlighted by the above mentioned 
rates. Last but not least, the error-prone rate of a gradient is presented 
and proved to differ from the gradient of a rate. 

As a consequence of these formulations, various tensor equations 
(e.g., equilibrium of elastic solids) can be differential (in time) as all 
tools are available. 

However, the validity of the tensor formulations of this section is 
restricted. As the time derivatives of the transformation tensor, Eqs. 
(2.3.11) and (2.3.12), are elaborated upon,* the existence of their inverse 
tensor, Eq. (2.2.24) is assumed. Hence, all cases where the inverse 
transformation tensor cannot be established are outside the scope of the 
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above mentioned rate operators. These are the cases where the 
dimensionality of the material coordinates is less than that of the spatial. 
For instance, unidiraensional (arches, beams) or bidimensional (shells, 
plates) material system of coordinates. Then the various curvature 
components ought to be considered. The conventional curvatures (namely - 
radii of tangent or normal circles) have to be accompanied by total 
curvatures, expressing the intrinsic properties. Rate operators for these 
values must be taken into account in an orderly fashion. 
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3. CONSTITUTIVE EQUATIONS 


3. 1 Introduction 

In Section 2, the equations, necessary for the precise treatment of 
constitutive equations, were presented. The present section contains a 
development of a three-dimensional theory which describes the incremental 
(rate) behavior of an elasto-thermo-viscoplastic continuum, in the presence 
of finite strain. Such a theory forms the necessary foundation for further 
developments of approximate theories for special structures (beams, plates, 
shells) . 

3.2 Incremental Theory of Plasticity 

The first formulations of incremental elasto-plasti c constitutive 
equations have been by given St w Venant (1870), Levy (1870), Von Mises 
(1913), Prandtl (1925) and Reuss (1930). These investigators have 
formulated various plasticity theories, in a style similar to that of the 
linear theory of elasticity, and ignored some important points, such as the 
correct stress rate and the difference between Lagrangenian coordinates and 
Eulerian coordinates in large deformations. A good historical review and 
additional details may be found in Refs. 53~54. It should be noted that 
these works form the basis of what might be called "an engineering theory 
of plasticity" and many works, books and researches rely upon them. 

The recent years brought a need for accurate solutions of 
elasto-plasti c problems with large deformations. It became clear that the 
classical elasto-plastic theories, mentioned above are not sufficient and 
it is necessary to give a more exact formulation of the constitutive 
equations. Three main research directions have been developed: The first 

direction begins with the works of Hill 1 7,55. i n this approach the elastic 
increment, which appears in the constitutive equations, is derived from a 
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finite elastic law. For that purpose a finite elastic energy function 
(potential) must be defined together with a finite elastic strain. 

A second research direction has been suggested by Sedov38 and 
independently by Green and Naghdi 56 . The elasto-plastic constitutive 
equations are based on thermodynamic considerations. The concepts of 
finite strain and a potential are also used in these works. The 
formulations are general and have not been applied to particular problems. 

The third research direction is not restricted to the theory of 
plasticity. Some thirty years ago new type of materials have been defined 
by Truesdell57-60 > the hypoelastic materials. The constitutive equations 
of these materials relate the stress rate with the strain rate. Some 
authors have shown that the hypoelastic definition includes, as special 
cases, various forms of plastic laws^) Detailed reviews on 

hypoelasticity are given in Refs. 39 and 65. 

The classical incremental theories of plasticity make use of an 
initial yield condition, a hardening rule, and a flow rule in 
characterizing the strain-hardening response of a material. Although these 
classical theories continue to be utilized extensively in finite element 
computer programs, this may be true only because more suitable models have 
not yet been developed. 

Comparison of the models with experimental results indicates 
relatively good agreement in uniaxial cases under simple loading 
conditions. However, for biaxial and triaxial cases and situations where 
the loading is cyclic, when creep and plasticity interact, and when the 
strain rates are high, the results are often in disagreement with 
experiment. 

The two most widely used yield conditions are the Tresca (maximum 
shear stress) and von Mises (J2 theory) conditions. For isotropic metals, 
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the von Mises yield condition generally provides a better description of 
initial yielding than does the Tresca condition. However, for rocks and 
soils, the Tresca condition is often used. Other yield conditions have 
been proposed, however, these have not found wide use because of their 
mathematical complexity. 

A flow rule is used to separate the total strain increment into 
elastic and plastic components. The most generally accepted flow rule, 
termed the normality condition, states that as the stress state of a 
material point comes into contact with and pierces the material's yield 
surface, the resulting plastic strain increment is along the outward normal 
to the yield surface at the point of penetration. 

The hardening rule provides a description of the changing size and 
shape of the subsequent yield surface, during plastic flow. In addition to 
simple expansion and/or translation, experimental evidence has shown that 
subsequent yield surfaces may exhibit corners, general distorsion, various 
Bauschinger effects, and dependence on prior cyclic history, strain rate, 
temperature and hold time, to mention only a few parameters^. For 
simplicity, most finite element programs make use of hardening rules which 
account only for expansion and/or translation of the yield surface. 

The classical isotropic hardening rule postulates that the yield 
surface expands uniformly during plastic deformation. In its simplest form 
wherein one assumes the von Mises yield condition and associated flow rule, 
the rate of strain hardening may be obtained by relating a value of 
equivalent total plastic strain to a point on a uniaxial stress-strain 
curve, so that a simple tensile test is all that is necessary to determine 
the hardening rule parameters. The simplicity of applying the isotropic 
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hardening rule has made it very popular in finite element plasticity 
analysis. 

In contrast, the kinematic hardening model of Prager-Z iegler 6 ? 
proposes that the yield surface translates as a rigid shape during plastic 
flow; the direction of translation being given by a vector connecting the 
current center of the yield surface and the current stress state. This 
gives rise to an ideal Bauschinger effect in which the reverse yield stress 
is lower by an amount of stress equal to the corresponding prior strain 
hardening. 

The Besseling-White (mechanical sublayer) model 68 makes use of a 
superposition of elastic— perfectly plastic stress states, in order to 
approximate strain hardening behavior. This model i3 often idealized 
mechanically as a parallel arrangement of elastic-perfectly plastic layers, 
whose yield stresses are adjusted to duplicate a piece— wise linearization 
of the uniaxial stress-strain curve (the number of layers being equal to 
the number of points selected on the stress-strain curve). Like the 
kinematic model, the mechanical sublayer model predicts a rigid translation 
of the yield surface. 

The hardening model proposed by Mroz 6 ^ employs the concept of a field 
of surfaces of constant work hardening moduli. Each point of a piece-wise 
linear uniaxial stress-strain curve is represented in stress space by a 
surface geometrically similar to the initial yield surface but of different 
size. The yield surface is assumed to expand and translate within this 
field, contacting and pushing each surface, along with it, as each is 
encountered. 

Krieg6g prop osed a two surface plasticity model, where the yield 
surface translates and expands within an enclosing "limit surface", which 
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also is allowed to translate and expand independently of the yield.- The 
hardening modulus is then assumed to be a function of the distance between 
the two surfaces at the loading point. Another model is the piecewise 
linear strain hardening theory of Hodge?®, which makes use of a yield 
polygon. 

As experimental evidence points out, isotropic and kinematic hardening 
tend to bracket actual material response in many cases and, for this 
reason, a number of combined isotropic-kinematic models have been proposed. 
Most models are based on a constant ratio of expansion to translation, 
although some results have been reported for a variable ratio based on 
accumulated plastic strain? 1 • 

3.3 Fundamental Assumptions 

In the treatment of elastic-plastic or elas t i c- vi s co pi as t i c 
deformations, we have to distinguish between the description as a 
thermo-mechanical process and the corresponding one by means of 
thermodynamic state equations. It is sometimes assumed that in the case of 
process which proceeds through non-equilibrium states, it is fundamentally 
necessary to start with a description of the pr o cess 1 9 » 3 7 » 7 2 . 
Alternatively, it has been proposed that one might assume local equilibrium 
for the elements of a body and therefore describe the state of the 
elements, in general, by state equations?3’‘75 . The consequences of 
adopting these two approaches become particularly clear when considering 
the influence of entropy. In the description of the process, entropy is a 
derived quantity and in- principle we can proceed without introducing it. 
In the description by state equations it is, on the contrary, a necessary 
state value, which, at least in principle, can be immediately determined. 
When restricting ourselves to homogeneous, quasi-statical thermomechanical 
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processes , the description by state equations can be viewed as equivalent 
to that by processes? 3 , 76 # The controversial issues will, thus, not be 
discussed further. 

We are dealing with large, non-isothermal deformation of solid 
polycrystalline bodies within the frame of classical continuum mechanics 
and thermodynamics. A phenomenological theory of such coupled 
thermo-mechani cal processes can be based on a material model including four 

different elements. 77, 78 

a) An elastic (or viscoelastic) element representing the reversible 
thermo-mechani cal processes governed by thermodynamical state 
equations; 

b) An element reflecting certain thermo-mechanical processes which 
lead to changes of the internal material structure, independent of 
plastic yield and thermally activated creep; 

c) A viscoplastic (plastic) element rendering also certain changes of 
the accompanying constrained equilibrium state; and 

d) An element representing thermally activated creep and relaxation 
phenomena, which also may be connected with corresponding changes 
of the internal material structure. 

Within the frame of the intended phenomenological theory, it is 
assumed that the respective thermodynamical state of each element is 
uniquely defined by the actual values of a finite set of external and 
internal state variables. Moreover, this enables one to introduce an 
accompanying fictitious state of constrained thermodynamical equilibrium, 
by means of a fictitious reversible process during which the internal 
variables are kept constant. This leads to a unique definition of 

reversible deformations77«78 > 
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Viscoplastic deformations are attributed to slip processes in certain 
crystal planes. These slip processes are based essentially on the notion 
of lattice defects. Roughly, we may distinguish local (bounded) processes 
restricted to the single grains (but occuring at the same time in all 
concerned grains) and global processes running through the whole body. The 
local processes comprehend the generalization and dissolution of lattice 
defects, the piling up of lattice defects at grain boundaries, etc. They 
are very sensitive to changes of the stress state, i.e. to the stress 
increments. The work Involved in these local processes is relatively small 
and mainly non-dissipative. The global processes are mainly dissipative 
and essentially governed by the actual stress state. The initiation and 
continuation of global processes, however, is always coupled with local 
processes . 

3. it Kinematic Considerations 

Basic to most of the postulated models of large elastic-plastic 
deformation behavior is the additive decomposition^ of d rg and E A0 (see 
Section 2.2.2) into elastic and plastic parts; 


d 


rs 







(3.4.1) 


The validity of this additive decomposition in the case of finite 
elastic-plastic strains has been questioned by -Lee and his 
associates 2 ^, 80-82 Lee’s 2 ^ approach is based on the total purely elastic 
unloading from the current state to an intermediate unstressed plastically 
deformed configuration, without any reverse or other kind of plastic flow. 
The major point in his theory is to decouple the total elastically induced 
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distortion and measure it from a relaxed unstressed state, which is only 
plastically deformed from the initial to the intermediate configuration. 
Accordintly, the deformation gradient F (our transformation tensor in 
Section 2) is decomposed in the multiplicative form, 

E P (342) 

F * F F ^ ' 

where F transforms a line element from the initial configuration to the 

E 

intermediate configuration, and F from the latter to the current 
configuration. The intermediate configuration is chosen in such a way so 

that F is unaffected by the presence of rigid body motion. The deformation 
rate tensors, d sr and d r3 are then defined. After some manipulations, Lee 
shows the following relation: 
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(3.4.3) 


where the subscript s denotes the symmetric parts. Generalization of Lee's 
theory for anisotropic elasticity was given by Mandel 8 ^. 

Lee's theory is based on the assumption that the elastic law does not 
change with the history of deformation and hence a total elastic unloading 
can take place. However, it has been shown 8 * 1 that after a fair amount of 
plastic flow has taken place, reverse plastic deformation will result soon 
upon unloading, even for small strains. Therefore, a total elastic 
unloading cannot have any physical significance. In view of this, the 
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theory of Lee appears as a special case of the theory of Green and Naghdi® 5 . 
Although not as general as the theory of Green and Naghdi, Lee’s theory, on 
the other hand, has the advantage of being more easily fitted with the 
physical property of invariance of elasticity with respect to plastic 
deformation. Mandel®3 in particular has pointed out that the Green-Naghdi 
theory is not convenient if one wants to include anisotropic elasticity 
effects. All this can be avoided by the use of the convected coordinates 
as proposed by Sedov3® and Lehmann®®. The formulation presented herein 
will follow the work of Lehmann. 

All quantities from here on will be related to the metric of the 
coordinate system u a in the deformed state. Hence, 
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( 3 . 4 . 4 ) 
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where the supercribed o relates to the initial state at time t. The 
deformation rate is. 




%{t 


-1 
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(f) 8 v - - yr') 0 . r 
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( 3 . 4 . 7 ) 


and the (♦) stands for the material derivative. 

The deformation gradient may be split into its elastic and its plastic 
components in the following manner: 
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(3.4.6) 


The use of capital greek subscripts and superscripts denotes 

parameters belonging to a fictitious intermediate state, defined by a 
fictitious reversible process with frozen internal variables, which is in 
general incompatible. The circumstance of the non-conti nuous configuration 
in the unstressed state has been observed by Sedov3&, who points out that 
convected coordinates, as used herein, become non-Euclidean in this 
configuration. 


fiditat rttonca >Mt 



Fig. 3.1 - Definition of Fictious Reference State 
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This multiplicative splitting of the metric change in the convected 
coordinates leads to an additive splitting of the deformation rate 
according to 


- syra */»{ Cf _ 1 (ff y } + sym ’/^(f" 1 )® (f)* r 


sym % Uf _1 )® r fj} - sym J /,{(f 1 )° (f"V g f*} 


(3.4.7) 


y y 

3.5 Elastic Deformations 

The present study is concerned with the structure of the constitutive 
relation of an elastic-viscoplastic (elastic-plastic) medium. The term 
elastic-viscoplastic means that the viscosity does not intervene in the 
elastic domain whose boundary, in particular, is well defined at every 
stage of the deformation. For simplicity we further assume that the 
thermo- elastic behavior of the body is isotropic and unaffected by 
inelastic deformation in the sense that the material 
constants characterizing the thermo-elastic behavior are Independent of 
inelastic deformation. 

It is convenient to work with the Klrchhoff stress tensor, t, in the 
current configuration, obtained from the Cauchy stress tensor by scaling 



(3.5.1) 
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where p denotes the current mass density, p Q the mass density in the 
initial state, and J the absolute determinant of the deformation gradient 
at the current configuration. 

The co-rotational stress rate also referred to as Jauman stress rate 
(see Section 2) will be 


a *a . .a > ."V a 

°B ■ ’.8 V d B "y 


a *a .a Y ,Y a 

r »t ■+ a t — a t 

g g y g g y 


(3.5-2) 


From Eq. (3.5*1), the following relations between the various rates of 
Kirchoff stress and Caushy stress are obtained 


♦a o *a , 7 .y a 
T “ — 0 + Jd P 

B P B .y B 


p o T .y a 

“ ff 6 + Jd .> a B 


(3.5.3) 


If a rate constitutive law is postulated between o ar >d d in finite 
inelasticity theories, then a potential does not exist, which is necessary 

in the variational or thermodynamics^based formulation of the problem. The 

y 

basic difficulty lies with the d.v-term. This is remedied by postulating a 
constitutive law between t and d. 

Thus, we can obtain a unique relation between the elastic deformations 
E 

represented by f® , the Kirchhoff stresses, t®, and the temperature, 
t 85,86. 
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( 3 . 5 . 4 ) 



Therefore, the thermo-elastic deformation is not independent of the 
inelastic deformation occurring at the same time. This follows from the 
fact that in the integrated form of the thermo-elastic stress-strain 
relations, Eq. (3.5.4), the stresses and the strains are referred to the 
deformed state of the body. 

In view of the present discussion and the discussion in the previous 
subsections, the hyper- elastic behavior described by Eqs. (3.5.4) and 
(3.5.5) will be replaced by a hypoelastic law. The hypoelastic law is a 
path-dependent material law, since it cannot be expressed in terms of an 
initial and a final state; it depends on the path connecting these states. 
Otherwise, if we did not make such a change, it would be necessary to 
retain the finite deformation measure in the constitutive law. For small 
elastic strains, there is practically no difference between hypoelastic and 
hyperelastic laws, as shown, for example, by Lehmann®®. 

The above could be illustrated by the following example. From the 
frequently used elastic stress-strain relation 
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we get 


A ■ 4c * -rf f ° 

which may be replaced by 

E y y 

.o 1 i a v & ,cn . £ ,a. x a 

d y - 15 K ■ 77^ t 8 V + aT S + aT 6 y 

We assume that inelastic deformation occurs if and only if 
F( , T, k, . . , Qy , • • . , A-yj * • • • ) “ 0 


(3.5.7) 


(3.5.8) 


and 


(3.5.9) 
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a 9F ♦ 

T v + ar T > 0 

a t 9T 
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for elastic-plastic material, 


or 


F(t®, T, k a® > 0 


(3.5.10) 

r ?o 

for an elastic-viscoplastic material. 

The function F represents the yield condition which bounds the domain 
of pure thermo-elastic behavior, in the ten-dimensional space of stress and 
temperature. The inequality, given by the second of Eqs. (3.5.9), is the 
loading condition. The actual form of the yield condition for a given 
material is determined by a set of so-called internal parameters, which are 
scalars and/or tensors of even order. The current values of the internal 
paramters depend on the initial state of the material and the history of 
the thermo-mechanical process. 

3.6 Thermodynamic Considerations 

Restricting ourselves to elementary processes, we need not analyze 
whether the applied heat arises from heat conduction or from heat sources. 
For the same reason it is not necessary, in our case, to introduce the 


61 



temperature gradient in addition to the temperature, and the body forces in 
addition to the stresses. 

The first law states, under our simplifying assumptions, that the rate 

of the specific internal energy, 0, is the sum of the rates of the specific 
mechanical work, W, and the specific applied heat, q: 


U - W + q 

The rate of mechanical work is given by 


( 3 . 6 . 1 ) 


W = — t° d Y (3.6.2) 

P o y a 

and may be split into an elastic and an inelastic part according to Eq. 
(3. 4. 7), 


w 
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W + W 


(3.6.3) 
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The rate of inelastic work must also be split into a part, W, which is 

S 

dissipated at once, and into another part W, which represents changes in 
the internal state. Thus, 


I . P v S D 

• 1 a T * * 

W • — t® d - W ♦ W 
p y a 
o 


Only W enters the entropy production 

D 

TS - q ♦ W 


The second law of thermodynami cs requires 


(3.6. H) 


(3.6.5) 


D 

W > 0 (3.6.6) 

We use as thermodynamic state variables the elastic strain, 


62 



E 

represented by f®, the absolute temperature T, and a number of other 

internal state variables (k, . . . , a®, . . . , A®^...), which may be scalars and 

E 

tensors of even order. The choice of f® and T as state variables is based 
on the fact that in pure thermo- elastic deformations, both quantities form 
a suitable set of thermodynamic state variables. The plastic strain and 
the total strain are unsuitable as state variables because, in general, 
they do not uniquely define the state of the material. A conflicting point 
of view has been expressed in Refs. 87-89. The remaining state variables 
are added for the sake of the description of the changes of the internal 
structure of the material. 

The specific free energy (Helmholtz function), $ , given by 

<|> ■ U - Ts (3.6.7) 

must be a unique function of the thermodynamic state variables 


E, 0 

“ <fr(f.y* T, k, a,y f ..., A^ . . . ) 


( 3 . 6 . 8 ) 


Since the elastic part of the deformation, according to our assumptions, 
does not depend on the plastic deformation, we may divide the free energy 
into two different components, as 


EE S 

<t> * , Tj ♦ $ ( T i 



(3.6.9) 


E 

where the first component, .refers to the elastic deformation and the 

S 

second, <t> , to the changes of the internal state. 

From Eqs . (3.6.1), (3-6.3), (3.6.H), (3.6.5) and (3.6.7) we derive 
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( 3 . 6 . 10 ) 


E S 

l = -sT + W + W 


Also, from (3.6.9), we obtain 
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( 3 . 6 . 11 ) 


By comparison of these last two equations, Eqs.(3.6.10) and (3.6.11), we 
may conclude that 
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( 3 . 6 . 12 ) 


For irreversible processes, this scheme of description has to be 
completed by some statements about the dependence of entropy production on 
the thermo-mechanical process. Under our assumption we need only deal with 
entropy production by dissipated mechanical work, in connection with 
inelastic deformation. Thus, we assume, in general 
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^ k > o 
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(3.6.13) 


where 
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(3.6.14) 


r aB r a $(r a t l- „ a a aB 1 

C T6 C Y5^ f 7’ T ’ • * a yt • • • a y6 J 

Eqs. (3.6.12) and (3.6.13) are the governing equations for non-isothermal , 

elastic-inelastic elementary processes. The specific free energy <t > , which 

determines the non-disspated work of the thermo- mechanical process, and 

the quantity , which governs the entropy production, must be specified 
according to the material behavior. 

3.7 Elastic-Plastic Model 

Elementary processes of elastic— plastic bodies may be considered as a 
sequence of equilibrium states, at least as long as the rate of deformation 
is moderate, so that the specific free energy is well defined in each state 
of the processes. We consider a simple example of isotropic behavior. In 
this case, the state of hardening can be described in the first 
approximation by a scalar state variable (beside the temperature). In 
order to simplify more, we assume here that the hardening is independent of 
temperature. Then the free energy can be written in the form: 


E 0 E E S , 

>(f“, T, h 2 } - <fc(f®,T) + *(h 2 ) 


(3.7.1) 


For the work during hardening It follows that: 



(3.7.2) 


W - <f*(h 2 ) - <J>(h 2 ) 

Provided now that there exists a uniquely defined relation (depending only 

D S 

on h^) between the dissipated energy W and the work during W - i.e. 
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(3-7.3) 
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then we have 
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[1 « o(h 2 )] *(h 2 ) - [t * c(h 2 )) *(h 2 ) 


(3.7.4) 


P S D 

The plastic work, W, as well as its parts W and W can, on the assumption 
made here, be represented as functions of the state variable h 2 . 

We can introduce a yield condition of the form, 

F(t°, k 2 ) - 0, (3.7.5) 

where is here a parameter characterizing the hardening. Since, on the 
other hand, the boundary of elastic behavior in the stress space can depend 
only on the state variable k^, then k^ is a function of h^, or 

p 

k 2 - k 2 (h 2 ) - k 2 (w) (3.7.6) 


holds. If the hardening is Isotropic, in a more restrictive meaning, the 
yield condition does not change its shape during plastic deformations, we 
have to put 

f[t®. k 2 ) - r(x®) - k 2 - 0 (3.7.7) 

Plastic deformation occurs if Eq. (3.7.7) holds and the loading condition 
is satisfied at the same time - i.e., 



o 


> 0 


We can use the state variable h 2 , for instance, by letting h 2 
Then it follows from Eq. (3.7.2) that 


(3.7.8) 

S 

- W. 
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(3.7.9) 


<j>(h 2 ) = h 2 
2 

Starting from the initial state h Q ■ 0, we get from Eq. ( 7 . 3 - -4 ) , 

P ? o 

W ■ [l + c(h )] h (3.7.10) 

P 

On the other hand, W can be represented from the yield condition, Eq. 

2 

(3.7.7), as a functin of k : 

P P ? 

W - W(k ] 

From Eqs. (3*7.10) and (3.7.11) we can determine: 

h 2 - h 2 (k 2 ) 

In the case of linear hardening with 

2 2 P 
k - k + 2B p W (3.7.13) 

o o 

we have 

o 

By putting c(h 2 ) - const. - c, thus assuming a constant relation between 
work of hardening and dissipated energy, it finally follows that 

2 2 

S , - k ~ K 

*(h 2 ) - h - — (3.7.15) 

2Bp o (1*c) 

For anistropic behavior and again, provided the plastic deformations 
are independent of temperature, the state equation for the free energy may 
be assumed approximately i.e. in the following form: 


(3.7.11) 

(3.7.12) 
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(3-7.16) 


<t> = <0 (f“, T) + *(h 2 , a“, 

3.8 Elastic-Vlscoplastic Model 

Thermo-mechanical processes in elastic- viscoplastic bodies cannot be 


considered as a sequence of equilibrium states, even in the case of the 


elementary processes, considered here. Elastic— viscoplastic deformations 

are associated with non-equilibrium states. One consequence of this fact 

is that we may get a continuation of a process without any change in the 

independent proces variables. This occurs, for example, in the case of 

creep with constant stress and temperature or in the case of an adiabatic 

stress relaxation under constant strain. In such cases, the body moves 

from a non-equilibrium state to an equilibrium state. 

In order to establish the constitutive relations for 

elas t i c- v i s co pi as t i c bodies, which in the limiting case becoming 

elastic-inviscidly plastic, we adopt the usual assumption that the 

stresses, which produce the inelastic deformation, may be expressed as the 

sum of the so-called athermal or inviscid stresses, , and the viscous 
* 

overstresses 
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( 3 . 8 . 1 ) 


This assumption, by no means, detracts from the "unified" concept. The 
rate- independent limit of visco-plastic constitutive relation was recently 
discussed by Travnicek and Kratochvil^ 0 . Hence, the total work rate can be 
partitioned in the following way: 
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The viscous part of the work is completely dissipated. Thus, we may write 

V D 

v 

* 9 

W = W (3.8.3) 


Regarding the plastic work, we have already stated that one part is used 
for changing the internal state and only the femaining part can be 
considered to be dissipated. Therefore, we must write 

PSD 

. . .P 

W = W ♦ W (3*8.4) 


So, we finally obtain 
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(3-8.5) 


w 

We have assumed that the changes of the internal state of the material 
can be regarded as a sequence of equilibrium states. Then, the specific 
energy is well - defined in each state of the process and we may take the 
usual overall statement concerning the specific free energy. In so doing, 


S 

however, we must be aware of the fact that into the part, W,. of the plastic 
P 

work rate, W, only the athermal stresses enter, since only these 
stresses are involved in the plastic mechanism. For the same reason, we 


can only introduce the athermal stresses into the statement concerning 


the dissipated plastic work 


9p 

w. 


On the other hand, we have to add the 


D v Dp 

dissipated viscous work W to W in order to obtain the total rate of 
dissipation. The different mechanisnms for determining the total 


dissipation and their coupling have been discussed by Perzyna^ 1 . 


We now consider an example in which the specific free energy has the 


following form 


69 


(3.8.5) 


EE S 

♦ - ♦ [ty. T) ♦ *(T t k, o“) 


E E 


*(f“. t) + k + f ( t) + ho“ a fl 


In this equation h denotes a constant with the dimension of a specific 
energy like the variable k and the function f(T). 

Furthermore, we assume that the dissipation is given by 
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where 5 < 1 and c denote constants. This leads to 


9 • • 9 p • 

w-w + W - (5 - 1 ) W - 5 ch a® + W 


Hence, we obtain 


( 3 . 8 . 6 ) 


(3.8.7) 


W ■ W - H • (i - c) W + 5 ch o° d Y 

t a 

On the other hand from Eq. (3.6.12) and (3.8.5) we have 
S 

' V 

W - k ♦ 2h o° a 
y a 

Eqs. (3.8.8) and (3*8.9) are compatible, for instance, if we put 

P 

k - (i - e) w 

and 
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( 3 . 8 . 8 ) 


(3.8.9) 


(3.8.10) 


( 3 . 8 . 11 ) 


From Eq. (3.8.10) it follows that, in our case, the plastic work, W, is 
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equivalent to the thermodynamic state variable k. This is still true if we 

take ^ as a function of k. But it does not hold in the general case when 5 
also depends on the other state variables T and a®. Eq. (3.8.11) shows that 

only in a very special case, a very unrealistic one, the state, variables 

are equivalent to the plastic deformation. 

From the thermodynamical considerations, it follows then that we may 

introduce the quantities k and a® , defined by Eqs. (3.8.10) and (3.8.11), 

P 

or any other equivalent set (W, cp o ho y ), as internal variables into the 
corresponding constitutive equations of the process description. 

The constitutive equations themselves are not yet determined 
completely by Eqs. (3.8.5), (3.8.6), (3.8.10) and (3.8.11). These 

equations only give the restrictive frame for the formulation of the 
constitutive equations. We may derive a complete set of constitutive 
equations, which is compatible with this frame, by the following further 
assumptions: 

a) the yield condition is of the form, 

F ■ (tj - c p o h a®) (t!J “ c P D ha °J “ “ 0 (3.8.12) 

- 

where t® denotes the deviator of the Kirchhoff stresses , 

b) the plastic deformation obeys the so-called normality rule. 
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(3.8.13) 


c) the relations between the viscous stresses and the Inelastic deformation 
rate are of the form, 


*<* ]_ !« , JL ft® - t®l 

d v “ “2^ S 2n l s V 


(3-8.1-4) 


d) the quantities £ an ^ c are constant. 


We can elimnate the athermal stresses x!J (which are not state 
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variables) from the equations of evolution by considering that the 
inelastic deformation can be expressed in two different ways. In one, the 
plastic mechanism is considered, and in the second, the viscous mechanism 
is considered. From Eq. (3.8.12) we then obtain, 

d* - 2i(t? - c p 0 ho°) (3.8.15) 

while from Eq. (3-8.14) we have 
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By comparing these equations for we get 
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(3-8.17) 


Following the course of the process in each state, the internal parameters 
p P 

W and a® and therefore also g 2 - g 2 (W, a y ) are known. Thus, we may 

calculate A from Eq. (3-8.17), and then all the other needed quantities 
P 

such as t^ and d^. 

This procedure fails at the point of transition from the elastic 


domain to elastic - viscopla3tic deformation domain, since in this instant 

P 

■ t® - T® and therefore A as well as d®, become zero. But in this case we 

P 

7 

may calculate d® from the following considerations: From Eq. (3-8.14) we 


obtain, because of d® - 0 and t® = t^ , 
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( 3 . 8 . 18 ) 
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On the other hand, we derive from Eqs. (3.8.14) and (3.8.15), observing 
that A - 0 and t® - t!J, 
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(3.8.19) 


From Eqs. (3.8.19), after multiplying by (t“ - cp Q h a o J, we obtain 
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( 3 . 8 . 20 ) 


Together with Eq. (3.8.18) this leads to 


a — hr ■ c p 0 h a V ' ^lt 


( 3 . 8 . 21 ) 
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Having the value of A, we may calculate d^ from Eq. (3.8.19), t y from Eq. 
(3.8.19), etc. 

3.9 Some Complementary Remarks 

Many thermodynamic considerations of non- isothermal, 
elastic-viscoplastic deformations refer essentially to the general 
fundamentals, which must be observed in describing such phenomena as 
thermo - mechani cal processes, and then discuss in particular which 
restrictions follow from the second law of thermodynamics. Only a few 


73 


papers attempt to describe completely such processes by state equations. 
Most of these papers introduce plastic strains, as thermodynamic state 
variables. But one may conclude from the consideration of the phenomena in 
the crystal lattice (dislocations, for example, which have completely 
passed through the crystal produce plastic strains but no changes of 
state), as well as from phenomenological observations (different states of 
hardening can belong to the same plastic strains), that plastic strains in 
general cannot be regarded as state variables. Furthermore, all these 
papers consider the plastic work as completely dissipated. This, however, 
is in contradiction with experimental results, from which it emerges that 
one part of plastic work is used for producing states of residual stresses 
in the lattice, which, when phenomenologically considered, cause hardening. 

The results, in the work presented here, can be extended to more 
complex constitutive equations by introducing more internal parameters or 
state variables, respectively. We may extend our approach to more general, 
anisotropic hardening materials by assuming (3ee Eq. (3.8.5)), for example, 
that 


EES 

4» - * (f°, T) + *(T, k, a®, A®J) 

- *<f“, T) . k • f(T) * 1% a’ a*. 


(3.9.1) 


Also, it may be more advantageous to replace the assumption in Eq. 
(3.8.13) for the plastic deformation rate by 


P V 

jO • 3F . n a& 6 


(3-9.2) 


This form of this model appears to be more suitable for representing some 
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experimental results in which second order effects and some deviations from 
the normality rule have been observed. Sometimes the normality rule is 
considered as a fundamental law based on an entropy production principle. 
But we should keep in mind that, since not all of the plastic work is 
dissipated, we cannot expect the total plastic deformation rate to obey the 
theory of plastic potential, even though the mentioned principles of 
entropy production are correct. 
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4. APPLICATION TO PROBLEMS IN EXTENSION AND SHEAR 
4.1 Introduction 

One of the most challenging aspects of finite strain formulations is 
to locate an analytical solutions with which to compare a proposed 
formulation. Typically, as a first problem, a large strain uniaxial test 
case is analyzed. The uniaxial tensile test is a common and simple way to 
characterize the stress-strain relation for a given material, since the 
tensor components used in the constitutive relation will have to be related 
to this uniaxial test. 

In Subsection 4.3 an example how the general constitutive relations 
developed in Section 3 can be applied to a particular material is shown. 
This material law is applied for all the other examples. 

The case considered in Subsection 4.4 examines the rate-dependent 
plastic response to a deformation history that includes segments of 
loading, unloading, and reloading, each occurring at verying strain 
amplitudes, for a bar. These are surely important problems to be 
considered; however, they only represent a partial test because the 
principal stretch directions remain constant. Finally, a problem which was 
discussed by Nagtegaal and de Jong 92 and others 93 as a problem which 
demonstrates limitations of the constitutive models in many finite strain 
formulations is the simple shear problem. This problem is solved as the 
last example. 
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4.2 - Uniaxial Irrotational Deformations 


The uniaxial tensile test is a common and simple way to characterize 
the stress-strain relation for a given material. Since the tensor 
components used in the constitutive relation will have to be related to 


this uniaxial test, and also to gain a physical understanding of the 
quantities involved in the analysis, it is both useful and instructive to 
express the tensor quantities previously discussed in terms of the uniaxial 
tension test variables. A homogeneous, uniaxial, irrotational deformation 
will be considered here as a first example. 

If the original length of the bar is i 0 and its present length is l, 
then, the transformation between the fix coordinate system and the material 
coordinate system is 



o 


and the metric is 



Ml 


“ g 


11 


11 



(4.2.2) 


(4.2.3) 



from Eq. (3.4.4) we find 



and 


(4.2.4) 


(4.2.5) 
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( 4 . 2 . 6 ) 


-1 il 


(f~M 


_o 

,2 


The component of the velocity vector in the fixed system are, 


1 


1 


l 


( 4 . 2 . 7 ) 


and in the convected material coordinate system they are; 

i i i i 

v * l V 1 “ .2 


( 4 . 2 . 8 ) 


The components of the deformation rate in the material system are 


then. 


l l 


.11 


i 2 i 

o 


11 


( 4 . 2 . 9 ) 


.1 


1 Jt 

It is worth while to notice that the material rate of a logarithmic strain 
is equal to the mixed components of the rate of deformation tensor. 

The unit normal vectors to the deformed and undeformed areas are one 
and the same unit vector directed along the bar axis, since the deformation 
is uniaxial and irrotational. Therefore, 
l 


1 o 
v " 1 


V 1 " l 


( 4 . 2 . 10 ) 


The force transmitted across the cross-sectional area of the bar is dP 

l 


dP 1 - — j- dP dP 1 - — dP 

o 


( 4 . 2 . 11 ) 


The corresponding traction vector components are: 

*1 P £ o P l 


VA 


1 


YA l 


( 4 . 2 . 12 ) 
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The components of the Cauchy stress tensor in the material system are 
obtained as: 



11 

o 


and the components 


_o P _£ P_ _o 

,2 A “ 


P A 
o o 


if P 

n 2 A 


P o \>£ 3 

o 


of the Kirchhoff stress tensor are: 


1 P p P _P £_ 

T 1 " p A “ A l 

K o o 


11 

T 


p p * 2 

o P o 
P a ,2 



(4.2.13) 


(4.2.14) 


T 


11 


p o P zf _P_ z£ 

p A * 2 ‘ A oZ 3 
o o 


Observe that the uniaxial component is the stress actually computed 

in most uniaxial tension test and is also inaccurately labeled some times 
as "true stress", since it is usually assumed to be equal to the "true 
stress" because p *■ p Q is satisfied almost identically for most metals in 
the plastic region. 

The Jaumann derivative of the components of the Cauchy stress tensor 
in the material system is: 
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d J f 1 | . 5 1 *t • 1 1 1 ,1 

at l'i J , 4 "i "i ' °i “i 


•1 • d fP i 

a * 1 dt A 


11 *11 A .1 11 11 .1 
; - a * d 1 c +o d^ 


*11 * O H 1 11 

a + 2 d.j a 


t 2 • l 2 Jl 2 

-fl- (l -2) ♦ 2 A (£ -2) - -2 -JL f£) 

dt l A £ 2 ; l x k % 2 ) g2 dt l A J 


11 


°11 ' d 1 a 11 


.1 * oj 

°11 d 1 ■ a 11 “ 2d 1 °11 


d fP t 2 p i fP l 2 1 if _d_ (-P-I 

dt l A £ 2 l A ^2 “ ^2 dt l A j 
o o o 


The corotationl (Jauman) rates of the Kirchhoff stress components can be 
similarly obtained: 


V 

T 


1 

1 




7 

T 


11 


_o _d_ fP 1_^ 

,2 dt 1 A i J 

1 o o 


V 

T 11 


if _d_ f_P 1> 

* ,2 dt 1 A 1 J 

1 o o 


o 

As previously noted in Subsection 2.5, the rate of the internal power 
can be expressed as a function of the Jauman derivative of Caushy stress 
and the rate of deformation tensor as in Eqs. (2.5.10) - (2.5.13). 
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4.3 ~ Example of Constitutive. Relations for an Isotropic Hardening Material 
For a carbon steel (45(DIN 1720) in a pure tension test at a moderate 
temperature and strain rate, we find the material behavior which is shown 
in Fig. 4.1^. 



Fig. 4.1 - Carbon Steel C 45 In Tension 
From this we may derive the stress-strain-temperature relations for loading 
in pure tension in the form 

o - o(f, T) (4.3.1) 

For our purpose it is more useful to write this relation in the form 
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p 

P c (T) W 

a - a(W, T) - - p + o q (T) 

c 2 (T) + W 

In our special case we get 

C 1 (T) - 72.42 - 36.03 • 10 _3 T , 

mm 

c 2 (T) - 7.35 - 8.04 • 10" 3 T , 

mm 

a Q (T) - 47.41 - 38.9 • 1(f 3 T , 

mrn 

with T in °K. 

We may consider the carbon steel approximately as an isotropic 
work-hardening material obeying the v. Misses-Hill yield condition. 
Furthermore, we assume that a constant ratio of 90% of the plastic work is 
dissipated. With these assumptions we get the following general 
process-description for the material under consideration in this example 
(see Section 3)» 

independent process variables: t®, T 

P 

P 2 P 

dependent process variables: (a) VI or k (W, T) f respectively 

(b) fj 


(4.3.2) 
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yield condition: 


F( 


V 


T, 


P 

W) 


t 7 - 
a 


2 ^ 

K (W ,T) 


0, 


K 2 (W,T) - | {o o (T) + 1 


P 

c, (T) W 


c 2 (t) + W 


loading condition: 


V V 

3F a . 3F • a _ 

— ; S 4 It t ' 2 S 

3 \ 


3k 

3T 


T > 0 


elastic strain rate: 

E 


"a lid v 6 t oi . ,o 

d y 2G 'S ” 1 + v t 6 6 y' 6 y 


(4.3.4) 


(4.3.5) 


(4.3.6) 


plastic strain rate: 

when Eqs. (4.3.4) and (4.3.5) are fulfilled: 


P 

d 


a 

y 



V 2 

2t e t 6 - ilsl 

6 3T 


3k 

P 

3W 


otherwise: 



(4.3.7) 


(4.3.8) 


rate of plastic work: 

P P v 

W - t® d 7 (4.3.9) 

7 O 

rate of applied heat (approximation; the exact formulation is given later 
in Eq. (4.3.17)): 
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p 

q « cT - 5 W 


(4.3.10) 


with 5 - 0.9 - const. c - 465 . ■ * „ (heat capacity) 

Kg K. 

In the process description the hardening parameter k^ depends uniquely 
P 

on the plastic work W and the temperature T. This fact leads to the 
following approach for the corresponding thermodynamics relations (see 

Subsection 3.6) to make sure that the plastic work is equivalent to 

thermodynamic state variables: 

free energy: 

EE S 

<t> “ <J>(f°,T) + h with h - * ; (4.3.11) 

entropy production: 

D ? y P 

W - 5(h) t° d - 5(h) W (4.3.12) 

t a 

From this approach and with connection with Eqs. (3.6.4) and (3.6.12), we 
get 


PSD P 

W - W + W « h + 5(h) W 

This means, as required. 



(4.3.13) 


P 

h - h(W) (4.3.14) 

In our special case, it holds that 5(h) - 0.9 - const. Therefore the 


S 

latent hardening « h becomes 

P P 

h - (1-5) W - 0.1 W 


(4.3.15) 
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From our approach, Eq . (4.3.11), we may also derive for this case, the 
exact relation between the rate of applied heat, the time derivative of 
temperature and the deformation rates. Using the Eqs. (3*6.5) and (3*6.12) 
we get 


D 

q + W 


e w - t s - t{- 


9S ..c* )3 • 3s *| 

— f v + 9T T + 3h hl 

9f° 


- t{- 


3 2 <t 


9f“ 9T 


V 2 2 

r a 9 d> i 3 d> 

V 9t 2 3h 9T 


h] 


(4.3*16) 


i.e. 


P 

q + E W 


2 E 2 E 
m [ 9 <& 9 

- T{— | T + E - y -~- 

9T 9f° 9 T 


V 

E 


In this equation 


(4.3*17) 


2 E 

- T i-| " T ) 

9T 


(4.3*18) 


represent the heat capacity at constant strain (d® - d“ - 0). We may 


consider this heat capacity as constant: 

E 

C(f“, T) - C. (4.3*19) 

Furthermore, we know from experimental results that the second term on 
the right-hand side of Eq. (4.3.17) can be neglected in most cases. So we 
can replace Eq. (4.3*17) by 
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(4.3.20) 


P 

q + 5 W » CT 

This is identical to the Eq. (4.3.10) which we used as an approximatin in 
the process description. 

4.4 - Uniaxial Cyclic Test 

The formulation described in Subsections 4.2 and 4.3 was applied to 

o 

the characterization of isothermic T - 400 k uniaxial cyclic response to 
loading program of variable strain amplitudes at a strain rate of 100 S~1 
(Fig. 4.2(a)). 

The dotted lines in Fig. 4.2(b) were obtained under the assumption 
that the total plastic work was dissipated. Due to the viscid effect, the 
rate-dependent stress-strain curves (solid lines) have continuous slopes at 
the shifting points. Note that the transient hardening causes the 
subcycles not to be closed. 
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4.5 - Simple Shear 


We take the same material as in subsection 4.3 and consider simple 
shear process as shown in Fig. 4.3. We denote this material as material A. 
The processes will 



Fig. 4.3 - Simple Shear 

be carried out, on the one hand, isothermically and, on the other hand, 
adiabatically . We find the solution of the problem by numerically 
integrating a system of first-order differential equations originating from 
the Eq. (4.3.4) - (4.3.10). In the first case (isothermic process), the 

total strain rate d° and the temperature T q are given, and in the second 
case (adiabatic process) the total strain rate and the vanishing of the 
applied heat are prescribed. 
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For comparison, we introduce, furthermore, a theoretical material 
whose yield condition is unaffected by temperature. This means, for this 
material, the hardening parameter k 2 is 

2 2 P 

k - k (w, T q ). 

In isothermic processes this material (denoted as material B) shows 
the same behavior as material A. But in adiabatic processes we have 
differences. For material A, the temperature Influences the yield 
condition as well as the elastic part of deformation. For material B, only 
the elastic components of the deformations are changed by temperature. 

Regarding this we must distinguish three cases: 

(I): isothermic processes with material A or B, 

(IIA) : adiabatic processes with material A 

(hardening rule depending on temperature), 

(IIB) : adiabatic processes with material B 

(hardening rule independent of temperature). 

The resuls for the shear stresses and the temperature are shown in Fig. 
4.i». We see that the differences between the shear stresses in the 
isothermic and in the adiabatic processes are mainly influenced by the 
dependence of the yield condition upon temperature. The differences betwen 
the cases (I) and (IIB) are negligible, but not the differences between IIA 
and IIB. It should be remarked that in the case (IIA), we get a maximum 
shear stress for v - 0.87. So for larger deformation we find in this case 
a softening effect due to the increasing temperature. With respect to the 
temperature the differences between the adiabatic cases, (IIA) and (IIB), 
are rather small, since the differences in the plastic work, in both these 
cases, are not so important. 
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The second-order effects are more influenced by the temperature than 
the first-order effects. This can be see from Fig. 4.5 The effects are 
partially' changed in the opposite direction (see stress o-j i ) . This is due 
to the strong influence of the temperature on the elastic deformations. We 
may conclude this from the fact, that the differences between the casses 
IIA and IIB are less than the differences between I to IIA or IIB, 
respectively. 

From other experiments, however, we know that the observable 
second-order effects cannot be explained by the influence of the elastic 
part of the deformation alone. We get more realistic result when we use 
stress-strain relations derived from the Eq. (3.4.2) with a small 
correction concerning the theory of plastic potential. In this case the 
difference between the second-order effects in isothermic and in adiabatic 
processes may be slightly less. But in any case, the influence of 
temperature on second— order effects is more important than the influence on 
first-order effects. 
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Fig. 4.4 - Response to Simple Shear Test 
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Fig. 4.5 - Response to Simple Shear Test 
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5. SHELL FORMULATION 


5.1 Introduction 

A basic framework for nonlinear or buckling analysis of thin shells by 
the finite element method is supplied by an incrementation of variational 
principles. Most of the nonlinear theories of shells are Lagrangian in 
character in that they employ as a reference configuration the undeformed 
state of the structure. In the construction of these theories it is common 
practice to start with a set of strain measures and strain-displacement 
relations (which are usually approximate in some sense), to introduce 
conjugate stress quantities, and to then postulate a variational principle 
of virtual work in terms of the variables of the theory. The next step is 
the incrementation of the principle of virtual work, which introduces 
further approximations and inconsistencies into the theory. 

A new procedure free from all the above cited limitations for the 
analysis of finite deformations, finite rotations, buckling and 
post-buckling behavior of arbitrary shaped shells of elastic or inelastic 
materials is presented in this Section. This formulation is a complete 
true abinito incremental theory. 

Subsection 5.2 describes the geometry and the kinematics of a surface. 
This formulation is used in the next Subsection which describes the 
geometry and the kinematics of the shell. The incremental principle of 
virtual power and incremental equations of equilibrium are introduced in 
Subsection 5. 


93 



5.2 Geometry and Kinematics of Surface 


The development of differential geometry of surfaces has been 

presented in numerous books , e.g,, McConnel £453 and Sokolnikoff [44], 
Following the notation of Ref. [44], a short and summarized development 

is presented herein. 

Surface s imbedded in three dimensional Euclidean space can be defined 
by: 

x 1 - x i [u a ] , i - 1,2,3 ; a - 1,2 (5.2.1) 

where x* are the coordinates of a point on the surface in the three 

dimensional system and u a are the coordinates on surface s. 

The covariant components of the metric tensor of the surface (or the 

first fundamental form of the surface) are: 


aB 


g x 1 x j 
g iJ a B 


( 5 . 2 . 2 ) 


The contravariant components, a^, of the metric tensor are defined by: 

a aB ' aB1 *1 <5 - 2 - 3> 

A unit vector normal to the surface can be defined as follows: 


1 aB J k 

n i " 2 E e ijk X a X B 


(5.2.4) 


_ O 

where £jj k and e are the permutation tensors for the space and for the 
surface . 


The total covariant derivative of a 
bring to the conclusion that: 


aB 


and use of the equality x 


i 

a; S 


x 


i 

B;a 
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0 


(5.2.5) 


g. . x 1 „ = 

B iJ a;B Y 


The implication here is there x 1 „ is orthogonal to x* and consequently to 
K a; B y 


the surface. Thereby, a tensor must exist so that: 

\-,e ■ b *e " l 


(5.2.6) 


Eq 3 . (5.2.6) are known as the Weingarten equations. With the help of 

tensor b one can build the second fundamental form of the surface. The 
aB 

total covariant derivative of the unit normal tensor is obtained by 
differentiation of the normality conditions: 


8 u \ " J 


0 , 


yielding 


i V.Y i 

n ■= - b x 

,a a y 


The integrability conditions on x 1 are: 


(5.2.7) 


( 5 . 2 . 8 ) 


a 2 X 1 
au° au 8 


a 2 x 1 
au 8 au“ 


(5.2.9) 


and they lead (from the covariant derivative of x „) to: 

a; B 


i _ i 6 i 

X a;6Y X a;Y6 “ .aBY X 6 ’ 


( 5 . 2 . 10 ) 


where R 5 is the Riemamn-Chr istoffel tensor of the surface given by: 
. aBY 


R 6 , aY - a6 + Y 6 - Y e Y 6 

• a8Y au 8 au Y aY e6 a8 eY 


(5.2.11 ) 
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( 5 . 2 . 12 ) 


R aBY6 “ a oe R .BY6 

Substitution of Eqs. (5.2.6) and (5.2.8) into Eq. (5.2.10) yields the 
Gauss-Codazzi equations : 


b aB;Y b aY;B 


(Codazzi) (5.2.13) 


b pB b aY b pY b aB R paBY 


(Gauss) 


(5.2.14) 


The Gaussian curvature of the surface k is defined by: 


k 


R 1 21 2 1 „ ctB Y6 

“a - 4 R aBY6 e e 


(5.2.15) 


where a - | a^l and e a ^ are permutation tensors of the surface (e 


aB aB 




or in alternate form: 
b 


K 


aB 


aB 


A b 

» a 


( 5 . 2 . 16 ) 


Another important invariant H, the mean curvature of the surface, .is 
defined by: 


H - 4 a° 6 b 0 (5.2.17) 

d aB 

Let us examine for awhile the role of the transformation tensor x* . For 

any surface contravariant tensor of first order t° one may write 

x 1 t° - T 1 (5.2.18) 

a 

Therefore, the components of t a are transformed to the space system. A 
similar transformation for the covariant tensor is by no means obvious 
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because the eq. (5.2.1) is not affine and there is no meaning to write u = 
u“(x i ). First, Let us define the inverse transformation tensor l“ by: 

( 5 . 2 . 19 ) 

This definition, Eq. (5.2.19), is obtained by multiplying Eq. (5.2.18) by 


L i ' 8 lj *6 a ° B 


g x| and by assuming that T 1 has only surface components. Then, 
i J 6 


. o T i 

T - t 

The inverse transformation tensor L® obeys the following rules: 

ij . a . B aB 

e 4 Lj . a 


( 5 . 2 . 20 ) 


( 5 . 2 . 21 ) 


L i X 6 


( 5 . 2 . 22 ) 


L® n 1 - 0 

Moreover, it can easily be verified that: 


(5.2.23) 


, a .a „ 

L i iS ■ “b "i 


<5. 2. 2*0 


» a e -- L ° n !s (5 - 2 - 25> 

The problems with the transformations, depicted by Eqs. (5.2.20) and 
(5.2.18), stems from the limitation on tensor T* to have only surface 
components. From a general tensor T 1 only the surface components are 
obtained by those transformations (in the case of higher rank tensor a 
little bit different component are lost). In order to overcome this 
problem, we define the surface system slightly different and add a third 
dimension by a unit vector normal to the surface. This system will be 
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called latter on T-system. 


The transformation tensors for the new system 


are defined as follows: 


and accordingly, 


'V 

a 


1 .2 

3 


(5.2.23) 


r L“ ; f - 1,2 


n i ; r - 3 


(5.2.24) 


The components of the co variant metric tensor for the new T-system are 


a . » v i „J _ _ *01 0 i 3 .3 

ta g ij r a aB 5 r 6 a i 6 r ^a 


(5.2.25) 


Similarly, the components of the contravariant metric tensor are: 


r a 


* ij 4 r L i 


.oB 


Now, the formal transformation of 
new t-system can be accomplished by: 


+ n i n 1 6^ (5.2.26) 

any tensor from the space system to the 




x 1 T a + 
a 


n*T 3 


or 


(5.2.27) 


T i “ L [ fc r " L i t a * n i fc 3 (5.2.28) 

This system, just defined, will be employed in the discussion of shell 
structures . 

Let us assume now that surface s is deformed in time so that the space 
coordinates of a point on the surface are time dependent. Hence instead of 
Eq. (5.2.1), we may now write 
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(5.2.29) 



According to the three dimensional formulation in section two , the x* 

system is the system at rest and the u a system is the material system (u a 
being independent of t). At this stage of the formulation we need to 
obtain the time derivatives of all the static parameters defined earlier in 
this section. 

The rate of change of the transformation tensor x^ is paramount in 
the description of the deformation of the surface. Hence repeating the 

results of section 2.3 [Eq. (2.3.12)]; 


Dx* 

a 


— (/>) . r q „ v k x s 
sk “ 


v" 4 v q 


(5.2.30) 


Dt „ a ' ' sk a ;a » 

where the velocity vector in the fixed system is given by Eq. (2.3.6), or 


i dx* 3x* • Dx* 

V " dt 3t Dt 


(5.2.31 ) 


However , the components of the spatial velocity vector can be rephrased as 
follows , 


v 1 - x* v Y ♦ n* N (5.2.32) 

From this the expression, is obtained through covariant 

differentiation 


V 1 - xi(v Y - b Y W) + V Y ♦ W ) (5.2.33) 

;o Y v ;ct a Ya ;a 

Comparison of Eqs. (5.2.33) and (5.2.27) is suggesting the possibility of 
simplifying the formulation by using the x - sy3tem. To accomplish this, 
there is a need to formulate an additional parameter, the material 
derivative of the unit vector normal to the surface, n*. By time 
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differentiation of g^j n i nJ = 1, one can come to the conclusion that the 

material time-derivative of n*,, is orthogonal to n* itself and located 

on a plane which is tangential to the surface. Hence, one can represent 


Dn i 

— r-r by a linear combination of x , or 
u t a 


£§r - C a x j 
Dt a 


(5.2. 3*0 


Differentiation of the orthogonality condition, g,. x n J = 0, and use of 

1J a 

Eqs. (5.2.30) and (5. 2. 3*0 yields the following expression. 


8 Ij V .o " J * 8 ij \ cT 4 ■ 0 


(5.2.35) 


Rearranging the terras of and substitution of Eq. (5.2.33) into Eq. (5.2. 3^ ) 


yield the following expression for 


- a aY [b V B + W ] x^, 
oB ;a J Y 


(5.2.36) 


In order to shorten the formulation we introduce two new parameters, 


f V T - b Y W - a Y6 <|» a Y6 [v« - b W] (5.2.37) 
,a * , a a 6a ; a 6a 


Na V * »,« ' 


\X * aY * J 


(5.2.38) 


Next, the different tensor will be denoted in the t-systera. The components 
of the velocity vector, Eq. (5.2.32), can be written as, 


V 1 - xj v r 

Eq (5.2.30), can be extended to: 


(5.2.39) 
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(5.2. 40) 


Dx 


_r 

Dt 


; r 


v 1 + 
. a r 


The second term on the right hand of Eq. 


Dn 1 3 
“Dt 6 r 

(5.2.40) 


of [Eq. (5.2.23)]. Note that by employing Eqs. 
(5.2.38); one can write 


stems from the definition 
(5.2.33), (5.2.37) and 


; r 


[4 * T _ * 


n 1 pj - [x‘ p T )s 


( 5 . 2 . h 1 ) 


It is useful for the forthcoming development to calculate the material 


derivatives of the inverse transformation tensor L® (or in the 
T-systera). By definition through Eq. (5.2.19), and in the x-system, Eq . 

(5.2.24), time differentiation yields 


^ X T 

Dt “ Dt Dt °3 

- [-L^ i> a y ♦ p°] & T a - [l* p Y ] & T 3 (5.2.42) 

Consideration of the structure of Eqs. (5.2.41) and (5.2.42) suggests the 
possibility of simplifying the notation considerably by introducing a new 
tensor i|>j, , in the x-system, in the following way: 


TA 


g x ‘ 

8 ij r Dt 


(5.2.43) 


A detailed derivation of Eq. (5.2.43), taking into consideration the 
defintions given by eqs. (5.2.37) and (5.2.38), shows that, 


< Ya ■ V? 6 A * (-%) S A * 


S A 8 * 0 • S A 
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Hence: 


*r 3 = " *3r " “ p a 6 r (5. 2 .^ 5 ) 

and 

v 33 - 0 (5.2.46) 

Now, we may write Eqs. (5.2.41) and (5.2.42) as follows, 


and 




i 




(5.2.47) 


~d ■ ' L f *'.a <5 - 2 - M) 

In view of these simple expressions, Eqs. (5.2.47) and (5.2.48), one 
can appreciate the importance of the introduction of the x-system. It is 
worth noticing, at this point, that, although most of the tensor quantities 
will be expressed in T-system, the mathemtical operations, like the 

material derivative, will be performed in the material system u a . This is 
done because the T-system is not a true material system, because of the 
third direction which is taken to be normal to the surface. 

The next step is to calculate the material derivative of the metric 
tensor in the T-system, Eq. (5.2.25): 

Da 

LA , - 2 - fg x 1 ! (5 2 ^9) 

Dt Dt Lg ij r A J Ko.d.W) 

Since gjj is independent of time, we have 


Da 


rA 


Dt 


«u [ 


DxJ, 

Dt 


Dx^ 

♦ x 1 ] 

a r Dt J 


(5.2.50) 


But from the definition of 
(5.2.25) we obtain: 


TA 


[Eq. (5.2.43)] and from Eqs. (5.2.47) and 
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Da 


_TA 

Dt 


l*fA * 


V 1 


(5.2.51 ) 


It is clear from its definition, Eqs . (5.2.43). that ^ is not a symmetric 
tensor. The description of the components of v rA in the u“ system is 


Da 


rA - (v „ ♦ v„ - 2Wb J 6® 
^ a; 8 8 » a a8 r l 


Dt 


(5.2.52) 


and the remaining of the (derivative) components, which permit r - 3 or A - 
3 or r - A - 3. vanish. 

The material derivative of the contravariant components of the metric 
tensor is obtained in a similar manner from Eqs. (5.2.26) and (5.2.48); 


^ - - a r V A ^ ) 


Dt 


.ft T ft 


- = ri> 


(5.2.53) 


- * V Y;p ’ 2 W 6 Yp> *« * 


r .a 
6 


The derivative of the metric tensor determinant a = det (a a5 ), in the u 


system is obtained by: 


Da A _3a 

Dt - at 


3a 


3a 


aB 


„ Da „ 
3a aB 


3a 


aB 


3t 


3a 


aB 


Dt 


(5.2.54) 


where it is recognized that a a& is a surface tensor so that its material 
derivative is f in fact, a partial derivative. 


— i 3 1 in fact , the cofactor of a given by: 
aB 
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3a 


da 


aB 

= a -a 


aB 


(5.2.55) 


Substitution to Eq. (5.2.54) and rearrangement of the components yields, 

(5.2.26) 


= 2a[v a - Wb“] - 2a /_ 
Dt ;a a J ,f 


Let us now introduce the rate of deformation tensor, d , by: 

r a 

d TA “ 2 ^TA + 


(5.2.57) 


and the spin tensor w by 

TA J 

“Va “ 1 Ka " ^Ar 3 ( 5 . 2 . 58 ) 

The rate of deformation and the spin tensors are the symmetric and 

antisymmetric parts of i^. The description of the components of d fA is 
simple, 


■fa - * »„« - 2 “». s ) *; 


(5.2.59) 


while the description of the components of u from Eq. (5.2.44) is 

TA 

“rA ■ 2 ! V aiS ' V B:J 4 r 6 a ' % 4 r S A * p s 4 r 6 a ■ 

- 5 ( Vb - V B,J 5 ? S A - < b ,T vr * ‘ 3 


* (b Br v> * “ iB ) 4 ? (5.2.60) 

Another important element in the description of the deformation of surfaces 


is the rate of change of the curvature, b 
[Eq. (5.2.6)] we obtain: 


aB 


From the definition of b 


afi 


b 


aB 



( 5 . 2 . 61 ) 
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The material derivative of Eq. (7.2.61) is then 


From Eqs . 
hand term o 


■Db „ Dn. . Dx . 

aB i i a;B 

s x + n — 

Dt Dt a;B i Dt 


(5.2.62) 


(5.2.6) and (5.2.36) it is readily obtained that the first right 
f Eq . (5.2.62) is zero. In order to calculate the material 


derivative of x* , there is a need to make use of the developed 

Ct J D 

connection between the time material derivatives and the total covariant 
derivatives (section 2.4). Following the derivation there and noticing that 
the space is Euclidean results into 


- V 1 (5.2.63) 

Dt ;aB 

Substitution into Eq. (5.2.62) of the covariant differentiation of Eq. 
(5.2.33) and then substitution into Eq. (5-2.13) yields the following 


expression for Dbg /Dt, 


Db 


-Si - nV 1 -b v Y ♦ b V Y + b V Y - b Y b W+W 
Dt i ;a8 Ya;8 YB ;a Ya ;B 8 Ya ;aB 


(5.2.64) 


or by using Eqs. (5.2.37) and (5.2.38); 


Db 


aB 


Dt 


b \i> + d 

8 V Ya p a;B 


(5.2.65) 


Another important variable, in obtaining material derivatives, is the 


gradient of the normal vector to the surface, n ^ . Once again from 
section 2.4 we have, for Euclidean space; 
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( 5 . 2 . 66 ) 


>g 

Dt 



;o 


From Eqs . (5.2.36) and (5.2.38) we obtain: 




+ 



(5.2.67) 


The final expressions in Eqs. (5.2.64) and (5.2.65) can be "translated", as 
they are, into the x-system (in opposite e.g. to Eq. (5.2.63)): 


Db 


fA 


Dt 


Db o 
aB *° 

— 6 r 6 a 


( 5 . 2 . 68 ) 
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5.3 Geometry and Kinematics of Shell 


5.3.1 - Geometry of the Shell Structure 

A shell is a three-dimensional body with one dimension small, compared 
to the other two. So, the shell structure can be viewed as some reference 
surface having material of some thickness on both sides. It is usually 
assumed that this reference surface is in the middle of the structure and 
so it is called "the middle surface." Hence, it is possible to describe a 
position of any material point of the shell by the following system of 
coordinates: 


f o 
u 


1 . 2 


(5.3.1) 


where u a are the material coordinates on the reference surface and and C 0 
is a material coordinate normal to the reference surface and assumed 

independent of time. It is possible to choose C Q as the distance of the 

material point from the refernce surface at time t - 0. Then the distance 

of any material point from the reference surface at any time will be 


constrained by: 

- h(u°,t) S C(e r .t) S + h (u°,t) (5.3.2) 


where h is half the thickness of shell in that particular point. Most 
generally, it is possible to describe the position of any material point in 
the shell by: 

x (e r ,t) - x l (u°,t) ♦ y l (e r ,t) (5.3.3) 

where x* is the position of the middle surface. The additional vector y 1 
can be decomposed, in the T-system, into a reference surface part and a 
part which is normal to the reference surface, or 
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(5.3. 1 *) 


y 1 (e r ,t) - x* y r = x* y a (e r ,t) + n i Y(e r ,t) 

The component of y*, which is normal to the reference surface, is the 
distance of the material point from the reference surface, meaning, 

e(e r ,t) * y(e r ,t) 

From Eq. (5.3.3) it Is clear that for e3 - - 0 we obtain: 

x(u°,0,t) - x i (u a > tj (5.3.5a) 

or, in fact: 

y 0, (u a ,0, t ] - 0 and Y(u a ,0,t) - 0 (5.3.5b) 

Without loss of generality, it can be assumed that, at t - 0, 

x i (e r ,0) - x l (u°,o) + Y(e r ,0)n i (5.3.6) . 

This means tht at time t - 0, the particular point is on the normal to the 
surface, so: 

y°(e r ,0) - 0 (5.3.7) 

Let us consider next a few approximating assumptions related to the 
character of the deformation. 

Assumption I : The material points which were on the normal before the 

deformation will be on the same normal to the surface after deformation. 

In other words: 

if:- y®(e r ,0) - 0 then: y°(e r ,t) - 0 (5.3.8) 

Substitution of this assumption, which is nothing else but the first 
Kirchhoff-Love assumption (see section 5.1), into Eq. (5.3*3) yields: 
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(5.3.9) 


x 1 (e r ,t) - x i (u°,t) + Yf e r , t )n i 

From Eq. (5.3.5b) it can be shown that it is possible to expand Y into a 
power series of £ 0 in the following manner: 


y( © r , t ) - l ii> n (u a ,t) 
n-1 




n-0 


(5.3.10) 


With regard to Eq. (5.3.10), it is possible to introduce yet another 
simplifying assumption. 

Assumption II : The shell is "sufficiently thin" in order to assume a 

linear dependence of x* on £ 0 then, 

Y(e r ,t) - E o * 1 (u°,t) (5.3.11) 

Substitution of Eq. (5.3.11) into Eq. (5.3.9) and use of the notation ^ f 
it yield: 

x i (0 r ,t) - x 1 (u°,t) ♦ E Q i|i(u 0 ,t}n 1 (5.3.12) 

It should be remarked that it is possible at this stage to assume a 
parabolic dependence on ? 0 and to get a theory which are similar to 
Reissner's theory. The linear dependence is the most popular in the 
different approximate shell theories (for this shells). 

Assumption III : The change with time of the gradient on the surface of the 

shell is negligibly small. A different formulation of this assumption can 
be found in Niordson’s work [95]. According to this assumption, points 
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which are on a surface which is parallel to the middle surface before 
deformatin will be on a parallel surface after deformation, meaning that 

is the independent of u a : 

SV.t) - x i (u°,t) + ^(tjn 1 (5.3.13) 

Additional assumptions lead us to the classical formulations. 

Assumption IV ; The change with time of a distance of a particular point 

I 

from the middle surface is negligibly small. This is in fact the third 
Kirchhoff-Love assumption and in our notation it means that 5 (or Y or it) 
are independent of time: 

x^e^t) - x i (u ot ,t) + n 1 (5.3.1**) 

Eq. (7.3.11*) is the common one (in the literature) for shell theories 
obeying the Kirchhoff-Love hypotheses. 

On the basis of the above four simplifying assumption, several 
formulations result, for the analysis of thin shells. These formulations 
are denoted below by capital letters. 

Formulation A : This formulation makes use of Assumptino I, only. 

Formulation B : This formulation employs Assumption I, and II. 
Formulation C : This formulation employs Assumption I, II & III. 
Formulation D : This is the classical thin shell theory based on the 

Kirchhoff-Love hypotheses or Assumptions I, II, III 
and IV. 

Formulation E : There exists a different formulation in the open 

literature 96 * Assumption I (accounts for transverse shear effects). This 
literature 96 * 9 ? which removes Assumption I (accounts for transverse shear 
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effects). This formulation starts by assuming that the shell is thin 


but y° * 0 in Eq. ,(5.3.4). In other words, they assume II without assuming 
I and obtain from Eqs. (5. 3.^) and (5.3.11). 

x i (e r ,t) - x i (u a ,t) + c o y i (u a ,t) 

- x i fu a t t) + E q [x* y a (u°,t) 4 'l>(u°,t )n i ] (5.3.15) 

Additional Formulation : It is possible also in Eq. (5.3*15) to impose 

additional assumptions about the shape of <l>, such as Assumptions III and IV. 
The kinematics of the shell structure can be developed on the basis of any 
of these assumptions (or others). The corresponding metric for the 
different formulations are described next, before developing the kinematic 
expressions. 

5.3.2 ~ The Shell Metric 

In this subsection we calculate the transformation tensors and the 
covarariant components of the metric tensor for every one of the shell 
formulations which were presented and discussed in the last subsection. 

First, for the case of the most general description (without any 
aproximating assumptions), from Eq. (5.3.3) we obtain the transformation 
tensor: 



Already in the notation of Eq. (5.3.16) there is a hidden assumption that 
the system at rest x* is cartezian. Otherwise, the transformation t*ensor, 

which in fact is 3x 1 /3e , would need an additional term which includes 
Christofel symbols of the space, on the right hand side of Eq. (5.3.16). 

From Eq. (5-3-^). one can write 
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Vi*i t»; - »i vi • »*[**/ . * k »; !0 * y !50 i «? 


i r w y 


a . [ i y 


(5.3.17) 

Similarly to the definitions given by eqs. (5.2.37) and (5.2.38) for v*-,r, 
we define here the following tensors, for y*;f 


a a _ a Y 

\s ■ y iB b s 1 


(5.3.18) 


y - b y + Y 
.B oB y ;8 


(5.3.19) 


By substitution of .Eqs. (5.3.17), (5.3.18) and (5.3.19) into Eq. (5.2.16) 
we obtain 

H\r - <4 K * 0 * ”\l 4 ? * IS y ;5 o * nl l! S ol s t (5 - 3 - 20) 


The covariant metric tensor of the shell is defined by: 

g rA * g ij *\r 


(5.3.21) 


Substitution of Eq. (5.3.20) into Eq. (5.3.21) and rearrangement of the 
terms yield the following expression: 

y 


* rt ■ l a o8 * l «s 4 **> * >. \ 6 * \ V *r 4 


a x 8 
A 


{ v + y + y Y } 6 a 6^ 
iy a;5o a y Y;£o a ;£o f °r °A 


(5.3.22) 


ly 8;50 * *B S:50 4 y B Xl 5ol 4 ? 4 J 


lv + Y Y I 5^ 5^ 

ty Y; 5 o y ;£o ; 5 o' °r °A 


112 



where a is the metric tensor of the surface, Eq. (5.2.2). The 
oB 

co n t r a va r i an t components of the metric tensor are defined from the 
following equality: 

l n • I A ° - 6 p (5.3.23) 

Similar to the surface formulation, it is possible to define here the 


~r 

inverse transformation tensor by: 


^i g ij ® xJ; A 


So that, 


rA -i , A 
L i x ;r - 5 r 


(5.3.24) 


(5.3.25) 


and 

i Ar - g iJ lJ Lj (5.3.26) 

It is worthwhile to notice that the covariant derivative of the shell’s 
metric tensor, Eq. (5.3.22), is not equal to zero. This means that raising 
and lowering of indices in the space of the shell is not permutable with 

r 

covariant derivatives. For example for any vector T : 


- (g 


r a 


tJ 

A ;< 


rA 

* g T 


A;* 


(5.3.27) 


Formulation A 

We shall formulate now the metric of the shell in view of assumption I 
(Formulation A). From Eq. (5.3.9) or by direct reduction of Eq. (5.3*20) 

for y° - 0, we obtain for the transformation tensor: 


-i 

x 


,r • 1-iK - ’0 * «? * I’.EO" 1 ! «? 


(5.3.28) 
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and for the metric tensor: 


g - la - 2b „ Y + Y 2 b„ b„ + Y Y „ } 6° 6 
B TA a8 aB Ya 8 ;a ;8 r 


o ,B 
A 


* iY ; » Y ;5 ot 6 ? S i 4 l Y ;S Y i 5 o> 5 ? 4 |l iC c Y ;5 ol S ? S » 

Formulation B 

When the second assumptin is imposed (Formulation B), we obtain from 
Eq. (5.3.12) for the transformation tensor: 

L ;r * { AK • 5 0 * b a ] + nl 5 o *;a J 6 r * { * 5 ? (5.3.30) 


r3 ,3 


-i 

x 


and for the metric tensor: 


k - la - 2E <l>b + E 2 Ui 2 b b Y + E 2 il» il> 16° 6^ 
g fA 1 oB o8 *o v Ya 8 o ;a ;B J r t 


U o * * ;a } + l?o * *. g } 6* 4 1^1 «J (5.3.3D 


3 t B , i.,2i ,3 r 3 


Formulation C 

When the third assumption is added (Formulation C) we obtain from Eq. 
(5.3.13) for x* : 

» J 

V " 6 r + { * ni} 6 r (5.3.32) 

and for the metric tensor: 

i rj ■ |a o 6 - 250 * b as 4 'o *\e> S T 4 <* 2 > s ? S A (5 - 3 - 33) 

Formulation D 

For the classical shell approximation theory, we obtain from Eq. 
(5. 3. U): 
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(5.3-3 1 *) 




5 b 
o a 


In 


*?♦ 


In 1 ! 




^ a o8 ~ 25 o b aB + Na "8 


,3 ,3 


b I> *r < ♦ «i 


(5.3.35) 


Formulation E 

Finally, for this formulation, we get slightly different expressions 
from Eq. (5.3.15). First, for the transformation tensor: 

x\ r - (xj [«I + Co J*] ♦ n*Co [9J} 6® ♦ {x* y y * n 1 *} 5^ (5.3.36) 

and then for the metric tensor: 




Co [ X 


+ 

a8 



*o l a 



y 

a 


V> s r 6 


B 

A 


+ { y a + *o [ K h + \ 6 t 6 a 
+ ^6 4 50 [ *B + V ]} 6 r 6 l 

♦ \y y y y * * z \ 6 ] (5.3.37) 

where: 

»° B ■ (?% - b°*) (5.3.38a) 

y - (b„ y Y * ♦ ) (5.3.38b) 

a >cr ;cr 


meaning that X and y are of the same structure as those defined by Eq. 
(5.3.18) and (5.3.19), but independent of Co. 
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5.3.3 ~ The Shell Kinematics 


The kinematics of shell can be fully described by the components of 


the rate of the transformation tensor , defined by Eq. (5. 2. M3), for the 

r a 

T-systera. From Eqs. (5.2.57) and (5.2.58) it is clear that defines 
exactly the strains and the rotations of the structure. If one is not 
satisfied with the internal description, there is a need to define the 
components of the velocity vector at every point of the structure. Similar 
to the metric analysis of the different possibilities to contract shell 
theories (see previous subsection) , in this subsection the components of 
the * tensor and the components of the velocity vector will be defined for 
every shell theory formulation of subsection 5.3.1. 

Following the kinematic analysis of surfaces, we define now the 
components of the velocity vector v by: 


-i dx* Dx* 
v " dt " Dt 

Substitution of the general definitions given by Eqs 
into Eqs. (5.3.34) yields: 


(5. 3. 34) 

(5.3.3) and (5.3.4) 


-i D r i i ci i v i 
V “ “Dt [x X a y n Y] 


(5.3.35) 


and through Eqs. (5.2.30), (5.2.31)* (5.2.36), (5.3*37) and (5.2.38) one 
can write 


-i i i a i D!‘ y i v i DY 
v - v * v Y * dT * P x 7 Y ♦ n ^ 


(5.3.36) 

where v 1 is the velocity vector of the reference surface. Eq. (5.3.36) can 
be shortened by making use of the T-systera as follows: 


-i 

v 


i i r i 0/ 

; r y r Dt 


(5.3.37) 
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We define now the rate of deformation tensor <J> for the general case Eq. 


-i 


(5.3.3). For this case the transformation tensor x is defined by Eq. 

5 * 

(5.3.20) and the metric tensor by Eq. (5.3.22). From the definition for 
Eq. (5. 2. ^3). we obtain 


’‘'r/i ’ g ij X ;F Dt 

Substitution of Eq. (5.3.20) into Eq. (5.3-38) and reordering of terms 
yield: 

y 


Dx^ 

-i _s A 


(5.3.38) 


T A 


f*-. + Ka *1 + *v 0 + Kx *1 ^ 


Di 


+ a 


6 , K 

— + a_. i — - 


a 0 a6 0 iB a Y6 a B aY Dt >6 0 Dt 


DY. 


* "Dt ‘ ’ "8 ’a * K >1 - \ 5 


* Kb y ,-E 0 * ♦« j Ie 0 * “M y is 0 * St y ,E 0 


W B r “ y B 


P 8 ) B y ’Eo " y B y ’Eo^ * P S I iEo^ 5 r S A * K> y % 


* *SY *6 y lE 0 * Vy * S« *o ~ St 2 * y o ~i Dt 2 


P X y ,E0 - 7 a y \o) - y ;Eo^ 6 T S ] 


+ ^Y6 y ;Eo y ;?0 + 3 Y6 


oy! 


j DY 

■tl 2 y i5 o ~ Dt 2 I 4 r { r (5 - 3 - 39) 


117 



where ii/ „ and p are given by Eqs. (5.2.37) and (5.2.38) and and 7 are 
aB a B a 

given by Eq. (5.3.38). Note that a is the metric tensor of the reference 

Op 

surface . 

For the sake of completeness we also include the expressions for the 


rate of deformation, d_, , and spin, w_., tensors according to the 

r is j I) 

definitions of Eqs (5.2.57) and (5.2.58), respectively, 

DA ft DA<5 

d fA " ^ d aB + d a« *B 4 d B5 A o + d 76 A a *B + 2 ^a6 ~Dt + *66 "St 


„ Dlf . Di 7 D7 D7 


y j 

fd v 7 ♦ d i 6 v 6 4 1 \a Py - ^q 4 a (-* B ■ y 

l a -yn O' . rr, fl vx *0 y.Cn 3 L a ov r\4- “vx i 


7B ;C 0 75 B ;£o 2 1 67 Dt 

D7. DY 


76 v Dt J ;£o 




^6 — Li2] + — £ y + y —i£2 l } 6® ■» 5 ^ 6^1 + 

B Dt ' Dt ;?o B Dt Jf 1 T A f A ; 


Dy Y 

{ H y^ + a L§2 y^ + Y 

' Q VX y.Crx y.r r> a vx nt J.r- - 


DY 


75 J ;£o J ;$o 76 Dt J i£o i£o Dt 


i£2i *3 ,3 , 
Dt 6 r 6 a 


, Di „ Di a Di „ Dr a D7 

- { d ♦ 1 f— 2i ♦ _i2_ 4 — 2_ j + ^ _Ii ♦ Y + 

1 aB 2 1 Dt Dt Dt oB a Dt a Dt 

+ y P> 0 II ,3 .8 , + % 4 j7 Py 7;l-0 

Y B Dt 6 r 6 A 2 Dt Dt y 7;Eo A B Dt 

. Y + Y fllESl (S 3 8 6 . S 8 J 3 ) . |!lliS2 ,■» . 

Dt ;£o B Dt 1 1 F A F A'' 1 Dt y ;5o 




(5. 3. ^0) 
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y y y 6 1 D> B DA a 

w rA = + W o> A 6 + W >6 + W >5 *a * B + 2 ~Dt " a 0Y "5t + 


DA fi A 

4 a -Y6^a “Dt “ “Dt + y 


dy dy 

— - - — - *y ] + y p - y p + 

a Dt Dt 8 J a ^8 8 


p « >6 * \ *a )! S T S l * Kb y ;'?o * %5 *8 y ; 'to * 


6 y 


!^. a [Af - »• fisj * , 

2 L a 8^ Dt a Y6 1 Dt y ;£o 0 Dt ; Dt ; £o 


D Y 

- y B -& 1 * P K2 Ko - y B y uo) * P B *1 - «! ft 


(5.3- 41 ) 


In Eqs . (5. 3. 40) and (5.3.41), d^ and w^ are the component of the rate 
of deformation and the spin tensors of the reference surface, defined by 
Eqs. (5.2.59) and (5.2.60). 

It is worth noticing here again that the expressions in Eqs. (5.3-39) - 
(5.3.41) are exact, but virtually impossible for someone to use them in 
numerical computations, because of their complexity. However, they may be 
used to check the accuracy and performance of the other formulations. 


Formulation A 

We develop now the kinematic expressions according to Assumption I 


*y 

(Formulation A), that is y * 0. The components of the velocity vector are 
obtained through differentiation of x 1 , Eq. (5.3.9), or by substitution of 

•y 

the assumption y - 0 into the exact expression of Eq. (5.3.36). 

Accordingly we obtain: 


-i 

v 


i 

V 


Y ♦ n -rr 


DY 

Dt 


(7.3.42) 
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ei 


The rate of transformation tensor is obtained also, by substitution of the 
right transformation tensor into an Assumption I [from Eq. (5.3.34) into 


the definition of Eq. (5.3*38), or .by reduction of Eq. (5.3.39)]: 


- j<t - it b 7 Y - ip b 7 Y + it b 7 b 4 Y 2 
fA 1 oB oB 8 Y8 a >5 a 8 


- f= _ a k< 5 vi v * 21) * v PY iB 

^ aB 5Y b a Y ' ' Dt Y b 8 Dt' Y ;a Dt 


- p Y + p Y -p Y T b 7 Y + b 7 Y 1 } 6° 6® 

w a ; 8 ^6 ;a 1 6 ;a a ;6 jr r A 


DY fl 

♦ { — L§. Y 

1 Dt 


;?o * b 8 Y Y ;?o + p 8 Y ;5 Q } 6 


3 ,8 
A 


(5.3.43) 


DY _ 

|y — i£2 

1 ;a Dt 


DY 

p b 7 Y Y - p Y } 6° 6^ + {Y ■ *5° } 6^ 6^ 
» -o ;?0 *o ;Eo' f A 1 ;£o Dt ’ T A 


The rate of deformation tensor for the Formulation A is 


Db 


d TA “ f d oB ~ Dt 


aB 


Y _ K 21+1 f 2 K Y 

Y b o8 Dt 2 Dt b Y8 Y 


Db 


2 . Pb YB W Y 2 * 


Dt o 


b’ Y j 


* b 7 b vo ♦ Y g ♦ 4 (X 


DY 


• 8 


DY 


4 y L2.1 1 jt a a b 

a Ore * Dt 2 l ";o Dt ;8 Dt Y 6 A 


DY DY 

+ 1 (_il y + Y f*3 ,8 + .8 .3^ 

.2 ‘TT Y ;Eo Y ; 8 ■ ' [6 T 4 A 4 f V 


DY .r« O o 

{Y + — ^2} 6^ 6 ^ 
1 ;£o Dt ' T A 


(5.3.44) 
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Similarly, the appropriate spin tensor is 


Db; 


Sa + Ke + Ko b 0 ‘ SB 0 Y + Ks b l b B + 5 (b 6o “Dt 


Db 51 DY .« 

- b.. )] * 2 ♦ i [Y - Y 


DY 


;o 


60 Dt 


2 1 ;o Dt ; 0 Dt 


p Y - p Y - p (bf Y Y - b 6 Y Y J} 6® 6? 
K 0 ;a a ;8 Y 0 ;a o ;0 T A 


DY DY 

+ {If — il y - Y — ) - P b 7 Y Y 
'2 l Dt 1 j Co ;0 Dt J ^Y 8 ;£o 


* "B Y .ect <5 - 3 ' ,,5) 

Study of Eq. (5. 3. 4*1) reveals that, in view of assumption I, there exist 
shear strains through the thickness. These strains are equal to zero in 
the reference surface [from Eq. (5,3.5b)], but different from zero anywhere 
else and their value depends on the distance from the reference surface. 
It is interesting to notice also that these components depends only on the 
gradients of Y and on the time rate of change of these gradients. 

Formulation B 

If we assume now a linear dependence of Y on £o, meaning that we 
impose the second assumption also, we obtain the components of the velocity 
vector through differentiation of Eq. (5.3.12), or by substitution of Eq. 
(5.3.11) into Eq. (5.3.42), for this case 


-i 

v 


y i 


D* i 


V - ^ * P x v 4 r TTT n 


’o ' ■- 'y ■» o Dt 
One should not confuse ij/ In Eq. (5.3.46), which is a scalar 


(5.3. 46 ) 
with the 


components of the rate of the transformation tensor, ^ . 
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In view of assumption II, the components of the rate of transformation 

tensor, i ., are obtain from substitutin of Eq. (5.3.11) into Eq. 
r a 

(5. 3. *»3). or 


Db 


-r r r y Y 3 > ^ D# 

— r U;hi) b + ii» b + a _ , J + b * tt 
V fA * aB %> l 1 aY B YB a aY Dt ' oB Dt 


Db: 


+ P B *;« ~ P a *;B^ + ^^YA b I b ft + b Y„ "*Dt^ 


Y6 a B Ya Dt J 


Y Dto 
+ b„_. b„ trr ♦ * 


Dip 


;B 


- p, *(bj *. o * t/ «; «j 


Ya B Dt ;a Dt 


Dt|> 


* l»p B - e 0 [» 2 » b p, - tt ' *)• 6 ? 5 « 


I - *P„ * «. I* 2 Pr * % sfll 5 ? ‘I * I* • 5?l 5 ? *2 <«•«> 


D'i'l 1 r o .3 


Wi x3 x3 


Similarly, the rate of deformation tensor, d rA , for assumption II, is 
obtain by substitution of Eq. (5.3.11) into Eq. (5-3- 44) , or 


5„ ■ K, ' * > oB §?! * #1 ** $ 


Dip 


Dip 


■>:s B *8-K-5i l **. B T? ,I K« 


a *B 
A 




t (5 - 3 


3 t B. ,e .3 


D* , ,3 ,3 


,H8) 


The components of the spin tensor u>^ f for this case, are obtained from 
subsitution of Eq. (5.3.H) into Eq. (5. 3.45), or 
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"r* ■ Kb * e o[*K 0 b l - % b a> 4 °e % ' p « *;B 1 


* •] »: b 8 4 ttb 


Db” Db Di|> 

6a ~Dt b 60 ~Dt ^ + V3 ^*;o Dt 


(5.3. ^9) 


Dl|/ f 
';B Dt 


*- -#) ■ <'{( b B * ■ b 0 * ♦ib^ 1 5 r s 5 


K 4 5 0 h 4 * - *„ f?> - » 2 s b p) ^ < - 5 ? ‘V 


In both Eqs . (5.3.48) and (5.3.47) it is customary to neglect the terms 


2 

which are multiplied by E , in accordance with the assumption that only 
the terms which are linear in 5 0 were kept before. Similar to the results 


under assumption I (Formulation A), here also the shear strains, d^j. , 
vanish on the reference surface. This strain components are called by 

Reissner [5.17] "the surface shear couples. The appearances of these shear 

couples causes surfaces which are parallel to the reference surface before 

the deformation not to remain parallel to it after deformation. This is so 

because surface gradient of does not vanish. It is customary to give a 

physical significance to the function t|> [5.58] through the definition: 


* - h ( u °> ll (5.3.50) 

h ( u° , 0 ) 

Clearly then, the term £ 0 i|> is the distance of a particular point from the 
reference surface, at time t. It is clear now, that the nonvanishing of 

the gradient of <l>, (* ) means that the change of thickness of the shell is 

inhomogeneous (nonconstant ) . 
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Formulation C 


If we superimpose now assumption III, that is, the independence of ^ 

of u° , then (as explained before) the surfaces which are parallel to the 
reference surface will remain parallel after deformation. This means that 

the change of thickness all over the shell will be the same. This does not 

limit the shell to be of the same thickness, throughout. According to 

assumption III, the corresponding velocity vector is obtained by 

differentiation of Eqs. (5.3*13). or 


-i 


v 



p 



, D* i 
5 o Dt n 


(5.3.51 ) 


This Eq. (7.3.51) is of course the same expression as Eq. (5.3.46) for 


Formulation B. In contrast though, th r ., which is obtained by a reduction 

i n 

of Eq. (5.3.47), has a much simpler expression. 


Db: 


J » {ip - p b T + ij> b 7 * a — -) 4 b I 

fA 1 c»B ^o LVl oV B >B o oV Dt' oB Dt J 


Db 




* I* "e * b B N 11 (5 r 5 a ' 5 f * f* 5? 1 6 t ‘a < 5 -3- 52 > 

Correspondingly, the components of the rate of the deformation tensor d , 

r a 

for this formulation are [from Eq. (5.3.48)3 


• K, - y * b oB If) * *5 tf *T? b T6 * ^ 0 


N 6 ‘ff)) { ? S A* )*lfl *?* 


D\Ji> 


(5.3.53) 
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Similarly, the components of the spin tensor w^, are obtained from Eq. 
(5.3.49): 


TA 


I w + r refill b„ “ b 11 
1 aB V K ya $ y& a' J 


* >1 b B * 


+ 


1 ° b 8 

l f h p _ 

2 l 6a Dt 


66 


Db 

-5? )»« 



6 

A 


+ i* Pg • S 0 [* 2 Py b 6^ ^ 6 r 6 a " 6 r ft (5.3.5V 

As expected, assumption III causes the surface shear couples to vanish. A 
formulation (development) based on Eqs. (5.3.52) ~ (5. 3. 5*0, with the term 
multiplyng S 2 neglected, can be found in Ref. 98. 

Formulation D 

The addition of Assumption IV bring us to the rate formulation of the 
classical shell theory, which includes all the Kirchhoff-Love hypotheses. 
Imposition of assumption IV, neglectng the change of thickness of shell 
with time (i|> - 1), into Eqs. (5.3.51) “ (5.3.54) yields the following 

expressions: 

for the velocity vector 


-i i _ y i 
v 1 - v - p x^ 


for the rate of transformation tensor: 


(5.3.55) 


Db, 


♦r, <\s ‘ 5 oi ( V« b e * % b I * VEt 11 * 


* <!«*»« b B * No “5t" 




a 


IP, - b P> (5 ? *! - 
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(5.3.56) 


and for the rate of deformation tensor: 


K 6 - N , * 


Db. 


IB K> 


*; * 


(5.3.57) 


TA 1 ~a6 '’o 1 Dt J ^ o l 2 1 Dt ">8 Dt 

the components of the spin tensor remain unchanged, and they are given by 
Eq. (5.3.54). 


A formulation with displacement vectors as given by Eqs. (5.3.55) - 

p 

(5.3.57), and with neglecting the terms which are multiplied by Eq. is 
the classical formulation and appear in many works (see subsection 5.1). 
Formulation E 

Let us formulate next the kinematics for Formulation E, for which the 
only requirement, from a physical point of view, is that the shell be thin 
and from a mathematical point that the expressions have a linear dependence 
on Eo* The components of the velocity vector are obtained through time 
differentiation of Eq. (5.3.15). 


-i 

v 


- V 1 + E (v 1 


-a A 
y 4 x 
;a a 


, Dv° Vi 
1 Dt ‘ p X 


i . ^ Di|/ i 'i 
V * * Dt " ) 


(5.3.58) 


The components of the rate of transformation, rate of deformation and spin 
tensors are obtained by an appropriate reduction of Eqs. (5.3.39) - 


(5.3.41), that is in every place instead of we substitute E 7° and 

p o 6 

instead of y we substitute Eo T , where and y are defined in Eq. 

a u a 0 a M 

(5.3.38). A kinematic development, based on this formulation and expressed 
in finite displacements and rotations, can be bound in Ref. So 
5.3.4 ~ Shift to the T-System 

The coordinate systems of the shell space x*, which were introduced in 
subsection 5.3-1. are not well adjusted for the introduction and 
development of the physical principles to follow. The reason for this is 
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that, on one hand the metric is confusing and of special character (like 
the nonvanishing of covariant derivative of the metric tensor) and on the 
other hand it is impossible to compare quantities which are defined in 
different coordinate systems. The solution to these difficulties is 
accomplished by shifting the different kinematic parameters to one 
coordinate system which is of a simple metric and obeys the Ricci Lema. 
Such a coordinate system is the T-systera which was defined in subsection 
5.2.1. 

The component of any space vector A* can be expressed in both systems, 
the T-system or the different x* systems, according to 

A 1 - xj a r - x 1 ;.^ (5.3.59) 

r A 

From this, we can obtain the connection between the a - components and the 
a - components, or 

a r - L[ x 1 ,a r (5.3.60) 

i ; a 

Following Ref. [98], we introduce the notation: 

vl L x 1 , (5. 3. 61a) 

A * l ; a 


-ir a r r i 
V A - L i A 


(5.3.61b) 


From their definition, Eq. (7.3.61), tensors p and p -1 are nothing else but 
euclidean shifters [99], and they will be used here to transfer the tensor 
components from the different shell systems to subsection 5.3.1 to the 
T-system. The shifters p and p -1 are absolute double tensors, and they 
obey the orthogonality conditions, 

— -in n 

(5.3.62) 


r -I a .r 
v a - v a “ 


r -1 £3 
v a ’ p r " 6 a 
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Also from Ref. [Q9], the metric tensor is transferred from one system to 
the other by: 


g 


fA 


7T <f> 

v r v a S* 


(5.3.63) 


in other words, the metric is conserved during the transformation. 

As mentioned earlier, the objective is to transfer the kinematic 
expressions of the last subsection to the T-system. A short study of 
subsection 5.3.3 reveals that the shift of the rate of transf ormation 


tensor, <», suffices to give all the needed internal information about the 
character of the deformation. The definition of 5/, from Eq. (5.2. A3), is: 



g ij 




(5.3.64) 


It is clear therefore that, in order to transfer i It is necessary to shift 
the transformation x and their material derivatives. We further note 
that we are not discussing a "legal" transformation from system to system, 


therefore it is impossible to perform a direct transformation of tensor 51 
without taking into account its components. 


From the definitions of x* and in the T-system, Eqs. (5.2.23) - 
(5.2.24) , and from Eq. (5.3.61) it is clear that the shift of the 
transformation tensor is taken place by: 



i A 
X A V T 


(5.3.65) 


For the material derivtive of x* , which is, in fact, the covariant deriva- 

J ^ 

tive of the velocity vector in the shell system, one must take a derivative 
of Eq. (5.3.65) as is, 


Dt ’ Dt (5.3.66) 
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From the chain rule and Bq. (5.2.47), one obtains 




i A , i Dv r 


Dt " V ;l\ v r + X /1 ~ Dt~ (5.3.67) 

Substitution of Eqs. (5.3.65) and (5.3.67) into Eq. (5.3.64) yields: 


*rA " g ij 


X 1 r v J u 0 + „i *' k a 1 

X ir V T V ; V T \ — J 


DP 


(5.3.68) 


Finally, use of Eqs. (5.2.43) and (5.2.25) results into the form 

a 


<4 


ra 


dp: 

V r W A ^irfi + a TT« v r Dt 


(5.3.69) 


The mixed components of J are defined by: 

* - g* r i r/( , 

and from the inverse operation, in Eq. (5.3.63). we have 

I* _ -H -IF T 

■ ■ a P . P _ <l> 


.A 


l TA 


(5.3.70) 


(5.3.71 ) 


Substitution of Eq. (5.3.69) into Eq. (5.3.71) and reordering of the terms 
with the help of Eq. (5.3.62) yield the following expression: 

n 0 

I* Q -H ,A -U Up A 

“ P A P J *.0 4 »» 0* DE (5.3.72) 


With the help of Eqs. (5.3.69) and (5.3.72) it is possible to express all 
the kinematic parameters in the T-system. The shape of the tensors v and 
its inverse, v ~' , will depend on the definition of the shell system, or, in 

fact, on the definition of x 1 , according to the different assumptions in 
subsection 5.3.1. 
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At this point, one may ask if there is a need for all the 
manipulations of the kinematics in the preceding subsections, since it 
suffices to find the expressions for y and y ^ and by these to uniquely 
define all the kinematics. However, inspite of the simplicity of the 
expressions in Eq. (5.3*61), it is not simple at all to find the 
expressions for the shift tensors, especially for y -1 . 

For the components of y we obtain a closed form expression by 


substitution of the general definition of x‘ from Eq.. (5.3*20) into Eq. 

9 * 

(5.3.61a), or 


v 


r 

A 





(5.3.73) 


where. 


Y A l 3 

a * .a 


3 A 

y - y 


(5*3*74) 


The expression in Eq. (5.3*73) is exact, and it is possible to introduce it 
into any one of the five formulations (A-E) or into additional formulations. 
As an example, we may substitute the classical assumptions (I-IV) into Eq. 
(5*3*73). In this case, y 1 - C 0 n 1 and 

i r - -V E 6 T ; y\ n - (5.3*75) 

.o a o a ;£o 3 

Substitution of Eq. (5.3*75) into Eq. (5.3*73) for y: 


vJ ■ f«I - 5. O 6“ 6 


a * r ♦ 6 3 *r 


o a ' V A A 3 


(5.3.76) 


which is the same expression as in Refs. [98] and [99]. 
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For the inverse shift tensor, p“1 , it is possible to make a similar 
development as for p, that is substitution of Eq. (5.3.20) into the 

definition of , Eq. (5.3.2H), and a repeat substitution into Eq. 
(5.3.61b) yields, 


p 



A<J» 




A rcJ 6 <t> + y r?co 


(5.3.77) 


In contrast to the expression for p, in Eq. (5.3-73). Eq. (5.3-77) Includes 
the metric tensor of the shell and is not expressible in terms of 
parameters for -r-system, only. In order to overcome this problem for the 
full Kirchhof-Love assumptions (Formulation E), Refs. [98] and [99] 
developed p _1 in terms of power series in Co, and kept only the linear 
term. If we assume linearity in Co, and kept meaning neglecting the term 
2 

multiplied by C Q in Eq. (5.3*35), and substitute Eq. (5.3.75) into Eq. 
(5.3.77), we obtain the same expression for the approximate p - ^ as in Refs. 

[98] and [99]. This expression is 


o‘ U r ' K * «• 4 


(5.3.78) 


At this state of our development, we shall not make any assumptions or 
approximations for the expression of p“i , and we shall keep the full 
expansion, as given by Eq. (5.3.77). 
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5.4 - The Principle of the Rate of Virtual Power 

5,4.1 - The Principle of Virtual Power and the Equation of Equilibrium 

The principle of virtual power (virtual velocities), as applied to 
shells, may be stated, in a similar manner as in the three dimensional 
f or mul at i ons , or 

| o U 5 ^ rA dV - | p B r 5v r dV - j y f 6v r - 0 (5.4.1) 

V V 

where all the parameters are referred to the material system of the shell 
space 

o rA - are the components of the Cauchy stress tensor 


* 


r a 


- are the components of the rate of transformation tensor 


§ r - are the components of body forces per unit mass 

F r - are the components of the external forces on the boundaries per 

unit area mass. 


v* - are the components of the velocity vector in the shell system. 

p - volume mass density 

y - mass density per unit area 

In order to simplify the first discussion, we limit ourselves to the 
case in which the external forces are applied only on the shell bounding 
surfaces (50 - ± h). Moreover, we assume that the shell is supported 
(constrained against translation) around the boundary of the reference 
surface in the transverse direction. A volume is defined in the shell 
system, Eq. (5.3.1) by 
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(5.4.2) 


pr i 2 ■? - 1 2 

dv - ig de de de g de de d? o 

where g is the metric determinant of the shell [{{? ■ detCg^)]. Eq. (5.4.2) 
can be written also in the following way: 

dV - da d£ Q (5.4.3) 

where da is a differential area at a distance 5 0 from the reference 
surface. 

Without losing generality, we can write for fS 

zu„) «.*.*> 

where g is the metric determinant in the x-system, and Z is a function of 
Cq , whose shape changes according to the definition of particular (r) 
system, or actually according to the different assumptions about the 
character of the deformation. In view of Eqs. (5.4.1) and (5.4.2), we can 
state for dV : 

dV - {F • ZUo) d ® 2 d e 0 - Z(!- o ) da d!- 0 (5.4.5) 

where da is a differential area on the references surface. In a similar 
manner we obtain the following expression for dA (differential area on the 
bounding surface) 


dAj - Z(± h) da (5.4.6) 

The expressions appearing as integrands in Eq. (5.4.1) need additional 
treatment also. It will be convenient to transfer all tensor quantities 
from the shell system to the x-system. For the first integral in Eq. 
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(5.4.1), by making use of Eqs. (5.3.61 ). (5.3.69), (5.3.62), by reordering 
the terms and by substituting Eq. (5.4.5), we obtain the following 

Dp° 

| o T * 6^ dV - j 6 “Dt 1 Z[ ^o ] d5 o da i5A ' 7) 

V V 

At this stage, we introduce the first approximatin by neglecting the 
2 

terms multiplied by ^ [inside the brackets of the Eq. (5.4.7)]. This 
linearization does not affect the term 6 ^ , which is independent of £ 0 . 

The approximation only affects the second term in the brackets, a term 

2 

which will always contain certain higher powers of 5 0 (but at least S q ). 

In view of this approximation, the term in the square brackets is at most 

linear in £ 0 . Hence, we can denote the coefficients of £ 0 by and the 

one which is free of £ 0 by $7^. In addition, we introduce the following 
expressions . 

+ h 

N** - | o r * Z U Q ) d? Q (5.4.8) 

-h 

+h 

M** - | • Z(5 o ) d? o (5.4.9) 

-h 

Note that both N 17 * and M 17 * are symmetric tensors, by definition. We can 
now write Eq. (5.4.7) in the following from: 

| o U Si n dV - | [N M t y tt . m’« « S ,]da (5.4.10) 

V a 

This form is used widely in the literature and it will be adopted, herein, 
in order to facilitate comparison with other works. 
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We next proceed to deal with the second volume integral in Eq. (7.4.1). 


Here also we transfer all the tensors quantities to the T-system. Vector F 


is transfered directly to the T-system in a standard way by employing the 


shift tensor y“ 
cons i der at i ons . 
express v* by: 


, but the transfer of 5Vj. needs some additional 
From the general definition for v*, Eq. (5.3.37), we can 




i 

A v 


(5.4.11) 


From this, it becomes clear that v are, in fact, the components of the 


velocity vector v in the T-system. The formal transfer of v^. to the 
T-system is therefore: 


Substitution of Eq. 


* 

Vj. - vjj Vj. (5.4.12) 

(5.4.12) into the second integral of Eq. (5.4.1), and 


(similar) transfer F 1 " into the T-systera, making use of the orthogonality 
relations of the shift tensors, Eq. (5-3-62), and of Eq. (5.4.5) yield the 
following expression: 

* 

| pB r $v r dV - j pB r 5v r . Z(e Q ) • dS Q da (5.4.13) 

V V 

TTA TT<fc 

We next define two new parameters, similar to N and M , by 

+h 

b r - | pB r • ZU Q ) • dS Q (5.4.14) 

-h 

+h 

d r - j pB r • S o • zU o ) . dS 0 (5.4.15) 

-h 

* 

Here, also we impose the assumption that in the expression of v, we can 
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neglect all the terms of E, 0 with powers bigger than one, whenever they 


ft 

exist. The term 6v r can also be separated into a linear component £ 0 with 
69^ as a coefficient and a components which is free of £ 0 , 60^°^. With 
the help of definitions given by Eqs. (5.4.14) and (5.4.15), the second 
integral in Eq. (5.4.1) can be expressed as 

pB 1 " 6Vj_ dV » | [b r 60^ + d r 60^] da (5.4.16) 

V a 

The treatment of the third integral in Eq. (5.4.1) is similar to the 
treatment of the second. Transfer of the tensorial quantities to the 
T-system and substitution of Eq. (5.4.6) yields: 


| y F r 6v r dA - j >F r 6v r [ZU Q ) 
A a 

We define now the following values: 


V h * 


] da 


f F - YF r [zU o ) 



] 

± h 


(5.4.17) 


(5.4.18) 





± h 


(5.4.19) 


By using these definitions and separation of $v as applied to the second 
integral, we obtain the following expression for the third integral: 


| 7F r Sv r dA - | [f r 60 ( ^ ♦ c F 60 ( J?] da (5.4.20) 

A a 
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Finally, from Eqs. (5.4.10), (5.4.16) and (5.4.20) it is possible to 
express the principle of Virtual Power, Eq. (5.4.1), in the following form: 

[ [N** &y + M** 6<- (b 7 + f F ) 69 ^ - (d F + c F ) ] da « 0 

J TT<p TT <p J i 

a 

(5.4.21 ) 

The equilibrium equations and the boundary conditions are obtaine from Eq. 
(5.4.21) by the application of the Gauss-Green theorem, where the character 

of the kinematics variables 6f x , 6k- . , 60^^ and 60^^ depends on the 

1T<p TTQ J J 

different approximations assumptions for the deformation of the shell 
structure. 

As an example we will derive next the equations of equilibrium and the 
boundary conditions for Formulation D in subsection 5.3.1; meaning, the 
classical theory based on the Kirchhoff-Love hypotheses. From the 
definition for ^ , in Eq. (5.2.44) and the expressions for the shift 

TT 

tensors for this case (Formulation D) in Eqs. (5.3.76) and (5.3.78) we 
obtain: 


sy 




<5lJ> 


1T<j> 



- v") *; 


6 


6 

* 


and 



k •; - •: 


5k* 




Db 


- jb 7 _ 6v + b 6v . 
a; 8 y a Y;B 


(5.4.22) 


- b] 6v v ♦ 6W „ + b, b„ 6W } 6° 6? 
8 a;Y ;o8 8 it 


(5.4.23) 


The expression for v in Eq. (5.4.11) reduces (for Formulation D) to: 
* 

3 


_ “ [v - £; (w 
T 1 a o ;< 


* \ \ 


)] a? 


♦ w6: 


(5.4.24) 


137 


From the above, one may write 


de^ - (dv ] 6® ♦ (6W) 6 3 r (5.4.25) 

r a j i 

d9 ( l ) - - ffiW + b 7 6v ] 6® (5.4.26) 

r ,o a y r 

Substitution of Eqs. (5.4.22), (5.4.23), (5.4.25) and (5.4.26) into the 
principle of Virtual Power, Eq. (5.4.21), making use of equations of the 
type: 

N® 6 6v . - (N oB 6v ) " N® & 6v (5.4.27a) 

a;B a ; 8 ;8 a 


,.aB . 

M b <$v. 


-O ~r„ ■ • 

and using the Gauss-Green theorem, for the first term on the right hand 
side of Eq. (5.4.27), yields the final expression of the principle of 
Virtual Power for Formulation D. 

} It- »“ 6 . 6 - b*). e . m’ ,6 . 8 b° - (b° . f°) . (d’ . o’) b“] 5» 0 




(5.4.27b) 


l -" 06 ** - " 08 So - "> 3 * f3 ) 4 * A,] 5 “l 3 


. | | [n« 6 - m’ s b° . h’ 0 * h’° b®] 


v„ • 6v » [M°® . d 6 * C 6 ] V. «« 

8 a 1 ;a J 6 


- M aB • v fl • 5W } ds - 0. (5.4.28) 

8 ;a' 

where v is a unit vector normal to surface s. Thus, the equilibrium 
equation for the classical theory will be: 
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- 0 • 


(5.4.29a) 


(N» e . b?)., - b^ p“ - m Y b» 

d“ B b * M° 6 „ . H oB b] b v 4 P 3 4 m Y - 0 (5.4.29) 

aB ;ap p t * 

P r - b r ♦ f r ; m r - d F ♦ C r (5.4.30) 

Moreover, the boundary conditions on s from Eq. (5*^*28), under the 
assumption that there are no forces applied on surfaces normal to the 
reference; are 

[n oB - b® + M^ a b B ] • Vg - 0 or 6 v q - 0 (5.4.31a) 


[M oB + m B l • v - 0 or 5W - 0 (5.4.31b) 

1 ;« •’0 

M oB » v » 0 or 6W - 0 (5.4.31c) 

f} ;a 

Similar to the definition of v, it is possible to define a unit vector t, 

tangent to s, and by this to break the partial derivative of <5.w into 

tangent to - s and normal to - s components 

<SW - SW • v a ♦ 6 W • t (5.4.32) 

;o ;y a ;s a 

Substitution of Eq. (5.4.32) into Eq. (5.4.28) and performing integration 

by parts, similar to Ref. [100], yield the final shape of the boundary 
conditions, Eq. (5.4.31). While Eq. (5.4.31a) does not change the next two 
do and instead of Eqs. (5.4.31 b,c) we have: 
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■ (m oB ♦ m B ) v 4 (M a6 v. t ) . - 0 or 6W = 0 (5.4. 33a) 

; a p p a ; s 

M aB • • Vg v q - 0 or 5W. a - 0 (5.4.33b) 

It must be noticed here that had we not assumed that there are no external 
forces on the surface normal to the reference surface then, the requirement 
for the vanishing of the static values in the boundary conditions, Eq. 
(5.4.31a) and (5.4.33). would be changed to the requirement that they would 
be equal to the external forces. 

Eqs. (5.4.29), which were obtained here for the classical theory, 
differ from the customary equations, as in Refs. [98] [100]. The reason 

for the difference lies in the definitions for N™ and M™. From Ref. 98. 
the stress resultant N and the stress couples M J are defined in the 
following manner. 

♦ h 

N™ - | • ZU Q ) d5 0 (5.4.34a) 

-h 


+h 

M™ - | a™ • Z o • Z(l- o ) de Q (5.4.34b) 

-h 

Moreover, similar definitions for the components of the body and external 
forces are given. Use of these definitions yilds the following system of 
equations [100]. 

‘ N ° 6 b af} * M aB b Y - M a8 ;oB - 0 (5.4.35a) 


- N aB „ - M YB b° „ - 


2(M YB b?) ;a - 0 


(5.4.35b) 


where Eqs. (5.4.35) were written without body or external forces. These 


equations, Eqs (5.4.35) replace Eqs. (5.4.29). 
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From the definitions, given by Eqs. ( 5 .- 4 . 3 1 *) > it is possible to 
developed equilibrium equations for the different formulations. For 
example, we obtain the equilibrium equations for assumption V of subsection 
5.3 by a similar process. Without external or body forces, the equilibrium 
equations for this case become 

- * » 33 N, * "° 6 5, - t»“ e *3,„e 


" 3S ^ b pY !6 * Ib B V * 21 6 ■ 0 


(5.4.36a) 


[-N° 6 + N 33 y 6 y a ♦ 


2M 76 Yb“] ;B - M 1 * Yb“ ;6 




_ { M >B y ° + lr u 3B „y .a A m 3Y 


y T y° ♦ M^ r 6 a]} 

2 L J J y y J 


|M 38 [y°(ifj ♦ ♦ y Y (ib“ . )“)])., - o 


(5.4.36b) 


N 0 * 3 - [M° 6 ♦ M 3S y°] ;B ♦ M 36 [Yb° ♦ ij] - 0 


(5.4.36c) 


N 33 Y - M oB b aB - [M 3$ Y] ;e ♦ M 38 [Y B - y y b u ) 


the related boundary conditions are: 

. B v # . (M°» Yv 6 t o ). 3 ■ 


or 6W - 0 


0 (5.4.36d) 


(5.4.37a) 


{2M >0 Y b“ + M 3> [y°(Yb5 + + ♦ A°)]} 


{[m 6Y y a + ~ y°(M 3e y Y + y 6 )] Vy t p } ;a - 0 or 6v Q - 0 


(5.4.37b) 
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V 6= ° 


(5.^.37c) 


[M° e + M 3tJ y°] 


36 


or 6 


u h 

Dt 


M 36 Yv 


6 


DY 

° r 6 “ 0 


0 or 6 ~ - 0 (5.4.37d) 

It will be also noticed that Eqs. (5.1*. 35a) and (5.4.36b) reduce for the 
equality to the known equations (5.4.35), as expected. Similarly, the 
boundary conditions, Eqs. (5.4.37a) and (5.4.37b), reduce to those found in 
Ref. [100]. 


5.4.2 - The Principle of the Rate of Virtual Power and the 
Rate of Equations of Equilibrium 

The principle of the rate of virtual power is obtained by total 
differentiation of the principle of virtual power [see Eq. (5.4.1)], 

Hit j si r, * - Hit j p B r dV - -jjj { IT F r «; r dA - o 

V V A 


(5.4.38) 

There is a need to perform a formal time differentiation of every one of 
the integrals in Eq. (5.4.38). Substitution of the definitions given by 
Eqs. (5.4.8) and (5.4.9), into the first integral in Eq. (5.4.38) yields, 

1 i Tt 1 { *r/i dv ‘ Hit j ["'H/i 4 "" SJ da 

V a 


or, 


I _ [ r dN™ + dM^ & 
1 j ^ dt 6y n dt . 6 


n , n™ . r1 FA *l**Tk\ , 

(f _ . + N 77 + M 77 — J da 


FA 


dt 


dt 


f f» r * S1 r » 4 M " 


d(da) 

dt 


(5.4.40) 
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where the time derivative of an elmental area of the surface is from Eq. 


(5.4.5): 

- -|jr (f? ) de 1 d9 2 = da (5.4.41) 

Substitution of Eq. (5.4.41) into Eq. (5.4.40) yields the following 
expression for I: 


I 


f 
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{[' 


fA 


dN 

dt 


JA 




6k 


rA 


+ 



dt 



<| 1 i ' c r< ) 

dt 


} da 


(5.4.42) 


The second integral in Eq. (5.4.38) can be developed in a similar manner. 
By using the definitions, given by Eqs. (5.4.13) and (5.4.14), we obtain 
[see Eq. (5.4.16)]: 


J, 4 4: ( f bF $e ( ? )+ d r da (5.4.43) 

i * at j j r 

a 

Finally, for the third integral in Eq. (5.4.38) we obtain from Eqs. 
(5.4.18) and (5.4.19): 

j 2 A _d_ | [ f f 69 <0> + c r 60 O)] da (5.4.44) 

a 

Use of Eq. (5. ft. 41) into Eqs. (5. 4.43) and (5.4.44), the addition of these 
two and use of definition, given by Eq. (5.4.30), finally yields: 


- l V * 59< r > 


r r d(M'" 

lP r ) * [» r l — sr 1 - I 


>0) 


dt 


dt 


(5.4.45) 



Hence, the principle of the rate of virtual power, generally, can be 
expressed by: 

I - J « 0 (5.4.46) 


where I express the rate of the internal stresses power and J the rate of 
the external (applied forces) power. 

Because of the nature of the elasto-visco-plasti c constitutive 
relations, there is a need to replace the total time derivatives in Eq. 
(5.4.42), by Jaumann derivatives. From Eq. (3.2.51) we obtain for I: 

. J.JA 


, f ifd r" r m $a . a M r* M rA i( fl -i 
1 " j ,l ~dT + N + N + N ♦.qJ 6 Va 


, r d J M rA . r u . a w r<j> . w r a ,o i r 

4 H — + v + M + M 6,f rA 


r r . dfair J 

[N FA ] FA 


dt 


n dffiif.J 

[«"] -ST- 1 - 


(5.4.47) 


As an example we shall develop now the principle of the rate of 
virtual power for the classical theory of shells. Formulation D in 
subection 5.3. For the sake of simplicity, we assume that there are no 
external forces and hence, the development deals only with the explicit 
expression for the internal power, I. It must be noticed that the 
addition of the integral which expresses, the influence of the external 
forces is not causing any difficulty at all. In order to obtain the 
explicit expression for I in Eq. (5.4.47), there is need to develop the 
expressions for the total time derivatives of 6Y_, and 6i c . Derivation 
according to time of Eqs. (5.4.22) and (5.4.23) yields: 
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(5.4.49) 


where we made use of the connections between the total time derivatives of 
covariant derivatives and the covariant derivatives of the total time one 
(see section 2). We define also the following symbol: 

A D7 


r Y ] : _°£ - f 

o Dt .a; 8 


(5.4.50) 

Substitution of Eqs . (5.4.22), (5.4.23), (5.4.48) and (5.4.49) into Eq. 
(5.4.48) and reordering of terms yield the explicit expression for I: 
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(5.4.51 ) 


As already mentioned, it is relatively easy to obtain the explicit 
expression for J, in Eq. (5.. 45), and thus complete the formulation of the 
principle of the rate of virtual power for the classical case of 
Formulation D). It must be noticed that the second integral In Eq. 
(5.4.51) is equal to the principle of virtual power itself and therefore 
can be erased from the expression for I. 

The principle of the rate of the virtual power is equivalent by 
definition to the equations of the rate of equilibrium and the 
corresponding boundary conditions. The equations of the rate of 
equilibrium can be developed therefore from eq. (5.4.51) by the same method 
we used to obtain the equation of equlibrium Eq. (5.4.29), from the 
principle of virtual power or by direct time differentiation of Eq. 
(5.4.29). In both case we obtain the following system of equations: 


r d J N o6 dV Y 0i d J M VB a op r a i H o0 r.Y i 

[—dr +— dTN J ;e ;$ N N l n J " 


dt t ,0 dt 


+ [ y ° 6 ] - >6 ^ 


^; 0-^0 


* ^ ^ - <6 M ° 6 3 ^;0 - M 


> »50 


147 







.ay 


* 


6 

• 0;5 


M 


75 




0 

.S;5 


- p°[n 3b 



0 


(5.4.52a) 
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2-2— b + 2JL- 4 2JL. b y b 4 N ° 6 {p 
dt a 0 dt ;o 0 dt 0 70 ,p a ;0 


* No *: B - 2 <o ^ * « 0 Vo s o 8 * 


2p .a b Y0 ~ d ct;Y0 + W .a;07^ M ;Y a d a0;5 


+ 2 M aB {** -dM- {M aB 

; 0 1 .Y;a a;Y ; 1 ;Y0 


O0 
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(5.4.52b) 


It might be thought that in the principle itself, Bq. (5.4.41), and in the 
equations of equilibrium, Eqs. (5.4.52) also, there exist additional terms 

which depend on N J °and M 3 . The requirement of symmetry of the Cauchy 

stress tensor in the system (-): 


TAir 


a ■ 0 , 


(5.4.53) 


can be transfer to the T-system in the following manner: 



Z 'U F rj> 




(5.4.54) 
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Integration of Eq. (5. 4.5*0 through the thickness of the shell and use of 
the definitions, given by Eqs (5.4.8) and (5.4.9), yield, among others, the 
following condition (for r » 3): 

f. v {N 37 ♦ M 3 V} - 0 (5.4.55) 
3>a o J 

From this it is clear that the terms in the brackets of Eq. (5.4.55) must 
be zero themselves without any connection to the permutation tensor. Study 
of Eqs. (5.4.51) and (5.4.52) shows that substitution of Eq.(5.4.55) will 

lead to the absence of N°^ and M a ^frora several terms (simplification). 


6. PLANAR CURVED BEAMS 


6.1 Introduction 

Section 6 deals with the kinematic and the constitutive equations used 
for the numerical analysis of planar curved beams. 

Formulation of problems concerned with finite deformation of beams has 
followed two different paths 1 O'*. Prescribing the beam by its deformed or 
undeformed centroidal axis and cross section, one may introduce at the 
outset beam stress resultants and their conjugate kinematic variables 
characterizing displacement and rotation of the cross section. Together 
with appropriate beam constitutive equations and a global balance law a 
consistent theory is obtained. Alternatively, one may imbed beam theory in 
the setting of deformable solid continua, in which, case one is concerned 
with local constitutive equations connecting the stress tensor with a 
strain tensor, which may in turn be expressed In terms of a combination of 
undetermined beam kinematic variables and functions of the beam coordinates. 
Momentum may then be balanced globally by integrating the local equations 
over the deformed beam configuration. Both paths will be considered in 
what follows. 

A complete abinito rate theory for the second path can be obtained by 
an appropriate plane stress approximation of the three-dimensional 
formulations presented in sections two and three. This approach is 
presented in Subsection 6.2. The rate form of the field equations for 
finite strains and rotations of curved beams according to the first path, 
can be obtained by a careful reduction of the shell theory, presented in 
section five. A simplified version of this formulation is presentld in 
Subsection 6 . 3 . Finally, three simple numerical examples are demonstrated 
in Subsection 6.4. 
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6.2 Plane Stress Approximation 

The constitutive relation presented in section three will be -rephrased 
here, in some different and compact form as follows 

p 

a) if F - (t* - C p Q g B*)(t* - C p Q g Bj) ‘ k 2 (W,T) = 0 (6.2.1) 

where, s,*, the Kirchhoff stress tensor, s* - — a*, and the temperature T 

9 k f K p K 

are independent process variables, t* being the deviator of the Kirchhoff 
i p i 

stress, s. , and W and (5. are internal parameters, then 
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(6.2.5) 

( 6 . 2 . 6 ) 
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( 6 . 2 . 8 ) 

(6.2.9) 

(6.2.10) 

( 6 . 2 . 11 ) 

( 6 . 2 . 12 ) 

( 6 . 2 . 13 ) 


P 

W - 0 (6.2.14) 

bJ - o (6.2.15) 

By definition, a body is said to be in the state of plane stress 
parallel to the u 1 , u 2 plane when the stress components a 1 3, a 2 3, a 33 
vanish^ 2 . It is well known in literature that the case of plane stress 
isdifficult to handle theoretically. Even linear elasticity has to 
treat this case in an approximate manner. To remove some of theoretical 
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difficulties Durban and Baruch 1 °3 introduced the notion of Generalized 
Plane Stress, where instead of dealing with the quantities themselves, one 
deals with their average values. 

In our case the problem is even more difficult. The nonlinearities, 
which the general three-dimensional theory takes into account will also 
cause a large change of the geometrical quantities in the u3 direction. 
Clearly, some assumptions are needed to treat the case of plane stress as 
a two-dimensional case. 

The first basic assumption is that the thickness, h, of the plate 
defined by the coordinates u 1 , u^ located in its middle plane, is small 
as compared with the other two dimensions. A second assumption is that the 
external forces act in the u 1 , u^ directions and are symmetrically 
distributed with respect to the middle plane. 

In a way similar to the procedure proposed by Durban and Baruch 1 
all the kinematic expressions are obtained by averaging the 
three-dimensional expressions. 

A basic assumption for the case of plane stress is that the components 
connected with the third direction are small and can be neglected. So, a 
new concept of generalized stress tensor is introduced 



It must be noted that in the linear theory of elasticity, where the 
geometry does not change, the averaged and generalized stress tensors 
coincide. 

So the three-dimensional incremental elasto-viscoplastic theory, 
developed previously, can be adopted for two— dimension plane stress 
problems. 
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6.3 A Simplified Version of Thin Curved Beam Element 



Fig. 6.1 - Reference Line of a Curved Beam 
A portion of the reference line for a curved beam is shown on Fig. 
6.1. The current arc length is denoted by s, while <fr is the current angle 
of inclination of the normar to the reference line, and p is the radius 
of curvature. 

The stress resultants acting on the beam cross section are the bending 
moment M, the axial force N, and the shear force Q. The external load, 
measured per unit of current length of the reference line, has the components 


p and p in the direction of the unit vectors e and e respectively. 

*s K n s n 

If v and v denote the velocity components in the direction of the 
s n 


unit vectors e and e respectively, the rate of extension is 
s n 

dv v 
d as p 


(6.3.1) 


The rate of rotation, <j>, of a given section is given by 


w » <t> 


3v v 

as p 


( 6 . 3 . 2 ) 
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while the generalized rate of deformation associated with bending is 

i 3v v 

r' • d<t> 0 ( n s > 

9s 9s 9s p 


(6.3.3) 


The rate of equilibrium equations for this simplified version may be 
put in the following form: 

|N_ Q |i. q|| + p + dp - 0 

9s 9s 9S H s F s 


fs + N £ + * H + Pn + d p n " 0 


(6.3. 4) 


fs + * + dQ * 0 


6.4 Numerical Examples 

The quasi-linear nature of the velocity equilibrium equations suggests 
the adoption of an incremental approach to numerical integration with 
respect to time. The availability of the field formulation provides 
assurance of the completeness of the incremental equations and allows the 
use of any convenient procedure for spatial integration over the domain B. 
In the present instance the choice has been made in favor of a simple first 
order expansion in time for the construction of incremental solutions from 
the results of finite element spatial integration of the governing 
equations. 

The procedure employed permits the rates of the field formulation to 
be interpreted as increments in the numerical solution. This is 
particularly convenient for the construction of incremental boundary 
condition histories. 

The capabilities of the presented models here-in have been evaluated 
through three simple numerical examples. The first example demonstrates 
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the capability of the plane stress approximation to predict deflections and 
stresses in a beam loaded by a constant moment. Figure 6.2 illustrates the 
beam and the finite element model. A quarter of the beam was divided into 
six elements in the vertical direction and into five elements in the 
horizontal direction. The external moment was introduced by six parallel 
forces acting on the section BC (see Fig. 6.2). 

The value of the external moment is 3500 kg/cm, and the material of 
the beam is CHR-17. The viscoplastic properties of the material were 
obtained experimentally from uniaxial tests in Ref. 9. This properties 
were collaborated into the present material model. 

The variation of the deflection of point E as a function of time is 
given in Fig. 6 . 3 . It is important to point out the value of the large 
deformation analysis. After ten minutes of the deformation is increased by 
42$ and at the same time there are important changes in the stress field 
(see Fig. 6.4). 
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The next example consists of a straight simply supported beam, loaded 
by a transverse concentrate force at the midspan. The beam is 25 inches 
long, two inches high and one inch wide. The material is stainless steel 
304 (Heat 9T2796). The material constants in sub section 6.2 were 
correlated with the uniaxial tension experimental results given in Ref.l04. 
The beam was subjected to a load of 2000 pounds at 1 1 00°F , this load was 
then held constant for 312 hr., and then increased to 2250 pounds at 
1 400°F . 

The primary purpose of this example is to compare the results, 
obtained by the two previously discussed models. The first one is the 
two-dimensional plane stress model, and the second one is the thin beam 
model as derived from thin shell theory. Figure 6.5 presents results in 
the form of load versus midspan deflection. The finite element model 
consists of five simple plane stress elements (dashed line in Fig. 6.5) or 
five sophisticated beam elements (full line in Fig. 6.5). 

It can be seen (Fig. 6.5) that the results agree quite well up to the 
312-hour hold period (points 3,4). During the hold period, the material 
hardens and only the beam model can represent this behavior after the load 
is further increased. 

The last example presents an analysis of a circular arch. The 
geometry of the shallow circular arch is shown on Fig. 6.6. The material 
is once again the 304 stainless steel. The arch is fixed at both ends and 
carries a concentrated load at the center. The el as to- visco pi as t i c 
analysis of this arch is performed with the aid of a ten curved beam elment 
model and with the inertia terms taken into account. The load P is assumed 
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The results of this 


to be applied In a quasi-static manner at t » 0. 
analysis are shown on Fig. 6.6, as the time-history of the midspan 
displacement. The response of the arch starts with the instantaneous 
elastic deformation at t - 0, followed by slow deformation up to point B, 
which can be considered as a limit point for the given value of the load P. 
Beyond point B, the displacements increase rapidly towards point C. This 
may suggest the existence of critical time for the prescribed load. 
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Fig. 6.4 - Stress Distribution 



Fig. 6.5 ~ Deflection vs Load 
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Fig. 6.6 - Circular Shallow Arch 



7. FUTURE RESEARCH 


As a consequence of these formulations, computational methods may be 
constructed. The incremental differential theorems lead to various finite 
difference methods. However, an integral theorem like the principle of the 
rate of virtual power calls for implementation by a finite element method. 
The discretization of the shell-like structure into finite elements and 
their systematic insertion into the integral theorems may yield a system of 
nodal motion differential equations. Numerous such applications are likely 
to be derived where large thermomechanical loads are anticipated. 

To develop geometrically nonlinear, doubly curved finite shell 
elements the basic equations of nonlinear shell theories have to be 
transferred into the finite element model. As these equations in general 
are written in tensor notation their implementation into the finite element 
matrix formulation requires considerable effort. The next effort will 
concentrate how to derive the nonlinear element matrices directly from the 
incrementally formulated nonlinear shell equations using a tensor-oriented 
procedure. This enables the numerical realization of all structural 
responses, e.g. the calculation of pre— and post— buckling branches in 
snap— through analysis and especially in bifurcation analysis, including the 
detection of critical points and the consideration of geometric 
imperfections. To avoid loss of accuracy care will be taken for a 
realistic computation of the geometric properties as well as of the external 
loads. Finally, the developed family of shell elements will be presented 
and its efficiency will be demonstrated by some applications. 
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